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ABSTRACT 

We probe the physical conditions in high redshift galaxies, specifically, the Damped Lyman-alpha 
Systems (DLAs) using neutral carbon (CI) fine structure lines and molecular hydrogen (H2). We 
report five new detections of C I and analyze the C I in an additional 2 DLAs with previously published 
data. We also present one new detection of H2 in a DLA. We present a new method of analysis 
that simultaneously constrains both the volume density and the temperature of the gas, as opposed 
to previous studies that a priori assumed a gas temperature. We use only the column density of 
C I measured in the fine structure states and the assumption of ionization equilibrium in order to 
constrain the physical conditions in the gas. We present a sample of 11 CI velocity components in 
6 DLAs and compare their properties to those derived by the global CII* technique. The resulting 
median values for this sample arc: <n(H I)> =69 cm" 3 , <T> = 50 K, and <log(P/k)> = 3.86 cm -3 
K , with standard deviations, <7 n (Hi) — 134 cm" 3 , <jt — 52 K, and o~i og ^p/k) = 3.68 cm" 3 K . This 



can be compared with the integrated median values for the same DLAs: <n(H I)> = 2.8 cm 
= 139 K, and <log(P/k)> = 2.57 cm" 3 K , with standard deviations cr n (# /) = 3.0 cm" 3 , <jt 



and a 1 



og(P/k) 



<T> 
43 K, 

= 0.22 cm" 3 K . Interestingly, the pressures measured in these high redshift Cl-clouds 
are similar to those found in the Milky Way. We conclude that the C I gas is tracing a higher-density, 
higher-pressure region, possibly indicative of post-shock gas or a photodissociation region on the edge 
of a molecular cloud. We speculate that these clouds may be direct probes of the precursor sites of 
star formation in normal galaxies at high redshift. 

Subject headings: Galaxies: Evolution, Galaxies: Intergalactic Medium, Galaxies: Quasars: Absorp- 
tion Lines 
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1. INTRODUCTION 

The high redshift neutral gas layers known as the 
Damped Lyman-a Systems (DLAs) are simultaneously 
well-understood and mysterious. On the one hand, the 
large SDSS survey has identified nearly 1,000 DLAs and 
produced a statistically significant description of the H 
I column density distribution function, the line density, 
and the contribution to the neutral gas mass density of 
the Universe, out to redshifts o f z ~ 4 (jProchaska et al.l 
[20051: iProchaska k Wolfe] [2009h . On the other hand, the 
precise nature of DLAs is still not well understood and 
basic physical properties such as the volume densities, 
temperatures, physical sizes, and masses are difficult to 
constrain due to the nature of their detection as absorp- 
tion imprints on background quasar spectra. 

As the dominant source of neutral gas in the Uni- 
verse between z=[0, 5], the DLAs represent a key link 
in the history of galaxy formation, as they likely pro- 
vide the source of neutral gas to fuel star formation at 

1 Some of the data presented herein were obtained at the W.M. 
Keck Observatory, which is operated as a scientific partnership 
among the California Institute of Technology, the University of 
California and the National Aeronautics and Space Administra- 
tion. The Observatory was made possible by the generous financial 
support of the W.M. Keck Foundation. 

Institute of Astronomy, University of Cambridge, Madingley 
Road, Cambridge, CB3 0HA, UK; raj@ast.cam.ac.uk 

3 Department of Physics, and Center for Astrophysics and Space 
Sciences, University of California, San Diego, 9500 Gilman Dr., La 
Jolla; CA 92093-0424 

4 Department of Astronomy and Astrophysics, UCO/Lick Ob- 
servatory; University of California, 1156 High Street, Santa Cruz, 
CA 95064 



high-z ()Wolfe. Gawiser. fc Prochaskal (|2003[ ). hereafter 
WGP03). The purpose of this paper is to constrain the 
physical conditions in DLAs by taking advantage of a 
relatively simple three-level system, the neutral carbon 
(C I) fine structure states. 

The ground electronic state of CI, 2s 2 2p 2 3 Pj, is split 
into three fine structure levels denoted as 3 Po, 3 Pi and 
3 P2- Following convention, we will refer to these fine 
structure states as CI, CI*, and CI**. The excited fine 
structure states are only 23.6 K and 62.4 K above the 
ground state, which makes them sensitive probes of con- 
ditions in cold gas. The advantage of using the fine struc- 
ture states of C I is that rather than determining only a 
line of sight column density, the relative excitation of 
the C I fine structure states allows for a determination 
of local physical conditions such as volume density, tem- 
perature and pressure. 

For almost two decades, the utility of C I fine structure 
transitions in probing cold gas has been recognized and 
applied to the high redshift Universe. The C I fine struc- 
ture levels are populated by collisional excitation and 
de-excitation, radiative decay, UV pumping and direct 
excitation by the Cosmic Microwave Background (CMB) 
radiation. Because of this sensitivity to the CMB, CI fine 
structure lines in DLAs were first used to determine the 
temperature of the C MB at high redsh i fts as a test of Big 
Bang Cosmology (i.e. ISongalia et al. I (|1994|) ). However, 
because the ionization potential (IP) of neutral carbon is 
below a Rydberg (IP C / = 11.3 eV < IP ff / = 13.6 eV), 
singly ionized carbon, C II, is the dominant state of car- 
bon in the ISM, and C I is not commonly found in DLAs. 
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Additionally, because the C I fine structure transitions 
are generally weak and suffer significant blending, it is 
only recently with modern high-resolution echelle spec- 
trographs that detailed measurements of the lines have 
been possible. 

C I has been detected in several high-z DLAs, often 
as a byproduct of a search f or molecular hydrogen (H 2 ) , 
i.e. iSrianand et al. I (|2005f ); lNoterdaeme et al. I (|2007bl ). 
Because they are photo-ionized and photo-dissocated, re- 
spectively, by photons of similar energy, the two species 
tend to be found together. Previous analyses of C I fine 
structure states in high-z DLAs have generally assumed 
a gas temperature in order to calculate the gas volume 
density and have generally found that for reasonable as- 
sumptions of the gas temperature, C I detections require 
relatively high densities. In this paper, we introduce 
a technique in which we make no a priori assumptions 
about the gas temperature. 

Inspired by the work of Denkins fe Shaval (|1979() and 
Uenkins fe Tripp I (|2001l ) who used C I fine structure ab- 
sorption to probe pressures in the local Milky Way ISM, 
we implement their technique o n high redshift DLAs. 
Most recently, Uenkins fe Tripp I ()2007f ) find that most 
of the CI in the Milky Way is in gas at pressures be- 
tween 3 < log (P/k) <4 cm' 3 K, with the distribution 
centered at P/k = 2700 cm -3 K. We will show that, in- 
terestingly, our r esults for the high reds hift DLAs are 
similar to those of lJenkins fe Trippl ()2007l ) for the Milky 
Way ISM. In addition, by invoking the sensible assump- 
tion of ionization equilibrium, we can use the population 
of the C I fine structure levels to constrain the total radi- 
ation field in some cases. This in turn acts as a check on 
the radiation field as determined by the 'CII* technique' 
(WPG03; WGP03) - a way of measuring the radiation 
field due to star formation in a DLA by equating the cool- 
ing rate measured via the CII* fine structure transition 
with the heating rate. 

This is the first in a series of two papers focused on de- 
riving the physical conditions in DLAs via the analysis of 
C I fine structure lines. In this first paper, we report the 
detection of C I in five high redshift DLAs, and analyze 
an additional 2 DLAs with previously published C I data. 
We derive limits on density, temperature, and pressure 
as well as the radiation field contained in the C I-bearing 
gas of each DLA. In the systems for which we have cov- 
erage of H2 , we compare the C I fine structure results 
to those derived by an analysis of the H2 rotational J 
level populations. In the second paper of this series we 
will extend our analysis to the C I-bearing DLAs already 
published in the literature, analyze the broader physical 
implications of the C I data and finally, propose a physi- 
cal picture of C I-bearing DLAs. 

This paper is organized as follows: We discuss our data 
and data analysis techniques in § [2j We present our pro- 
cedure for analyzing the CI data through the specific 
example of DLA 0812+32 at z abs = 2.626 in § 02 We 
then summarize the results for each of the C I-bearing 
DLAs in § [fj § [5] contains discussion of the DLAs pre- 
sented here and we conclude in § [B] Throughout this 
paper we make the standard assumption that the ratio 
of column density of element X is equal to the ratio of 

volume density, = ' we are analyzing 

average conditions along the line of sight. 



2. DATA AND METHODOLOGY 

The C I sample presented here represent serendipitous 
discoveries made during the course of an on-going cam- 
paign to obtain high resolution spectra of DLA targets 
for detailed study. Data for this paper was taken primar- 
ily with the HIRES spectrograph (Vogt et al. 1994) on 
the Keck I telescope, with a typical decker that resulted 
in an instrumental resolution FWHM = 6.25 km s _1 
. Details of specific observations are given in Table [TJ 
Data were reduced and continuum fit using the standard 
XIDL 5 packages. Because of the complex blending of 
many close fine structure transitions that typically ex- 
hibit several velocity components, we were not able to 
successfully apply the app arent optical depth technique 
(|Savage fe Semb ach 1991) for a measurement of the col- 
umn densities. We also did not hav e the high signal 
to noi se ratio and resolution of the Uenkins feT ripp 
(|2001[) study of interstellar C I, in which they develope d 
a modified AODM technique (j Jenkins fe Tripp I l200ll) . 
Instead we used the VPFIT package, version 9.5 6 to 
measure the CI fine structure column densities, N(CI), 
N(CI*), and N(CI**), redshifts, and Doppler parame- 
ters (b), where b~\/2a, and a is the velocity dispersion 
in km s" 1 . W avelengths and f-values were taken from 
IMortonl ([2003T ) 7 . 

In general, we simultaneously fit as many CI multi- 
plets as possible. In all cases, we included all multiplets 
that fell redward of the Lya forest, which generally in- 
cluded the strongest multiplets at 1656A and 1560A , and 
usually the multiplet at 1328A. In cases where a multi- 
plet fell within the Lyman a forest, yet did not contain 
any obvious blending with a forest line, we included that 
multiplet in our analysis. We rejected sections of spectra 
that contained obvious blending. 

To determine upper limits for non-detected fine struc- 
ture states, generally for N(C I**), we used VPFIT, which 
allows for an estimation of the ler upper limit by in- 
putting a "reasonable guess" of the column density and 
re-running the fit to give linear errors. The error is then 
taken as the la upper limit. For a conservative reason- 
able guess we used the measured column density of the 
line with the smallest measured column density, usually 
N(CI*), since it is almost always the case that N(CI**) 
< N(CI*). This results in a conservative upper limit to 
the N(CF*) value. 

We present one new detection of H2 in DLA 2340—00, 
and we analyze the H 2 in DLA 0812+32. While it is likely 
that the other systems contain measurable H 2 as well, we 
lack the spectral coverage to confirm this speculation. 

All other information about each DLA, including 
metallicity, dust-to-gas ratio, and log N(C II) was de- 
termined using the Apparent Opti cal Depth Method, 
AODM (jSavage fe Sembachl Il991h . unless otherwise 
stated (i.e. in the cases where we performed a full com- 

5 http:/ /www. ucolick.org/~xavier/IDL/ 

6 http://www.ast.cam.ac.uk/~rfc/vpfit.html 

7 Note, however, that there is some confusion over the correct 
f-values for the C I transitions given a more recent set determined 
by E. B. Jenkins (2006, private communication). Since we cannot 
det ermine which s et of f-values is more correct, in this paper we use 
the Morton ( 2003) values to be consistent with what has previously 
been used in the literature. 
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ponent by component analysis). To measure the log-^r^- 
we must employ the conventional method of measuring 
N(CII) by proxy using N(Si II), because the CII A1334 
transition i s saturated in every DLA presented in this pa- 
per. As in IWolfe et all (l200l we let [C/H] 8 = [Si/H] 
+ [Fe/Si]mt where [Fe/Si]i„t = —0.2 for a minimal de- 
pletion model or |Fe/Si1j n t = 0-0 fo r a maximal depletion 
model. While iWolfe et all (|2004D estimate the error m 
the measurement of ^rj- to be 0.1 dex, we use the more 
conservative estimate of 0.2 dex. 

Details of the observations are given in Table [TJ while 
the details of the C I measurements are summarized in 
Table [2] In the following we provide a brief summary of 
the data analysis of each DLA in our sample: 

•DLA 0812 + 32, z abs = 2.62633 : DLA 0812+32 
at z aos =2.6263 3 has been studied extensively (see 
IProchaska et al. 1 (|2003D ). and boasts one of the high- 
est k nown DLA metal li cities and H I column densi- 
ties. IProchaska et al. I (|2003| ) used the many available 
transitions to show that the relative elemental abun- 
dance pattern i s sim ilar to that of the Milky Way. 
IJorgenson et al.l ()200 9 ! ) presented the first direct evidence 
of cold (Tthermai < 78 K, (115 K, 2a)) gas at high red- 
shift using a curve of growth analysis of a sub-resolution, 
narrow velocity C I com ponent in this DLA, and recently, 
iTumlinson et al.l (|2010l ) presented the detection of HD. 

In this work we simultaneously fit data from three Keck 
HIRES runs with details given in Table [U Eight multi- 
plets were used to fit the C I fine structure lines; A1656, 
A1560, A1328, A1280, A1279, A1277, A1276, and A1270. 
The A1276 - A1280 multiplets were blended with a C IV 
doublet at z aos = 1.992. Figure Q] contains the spectral 
data of DLA 0812+32, in black, overlaid with our fit, in 
red, and the C I fine structure lines marked. Each veloc- 
ity component is denoted by a different linestyle (solid, 
dashed, dotted, etc.) while each fine structure state is de- 
noted by a different color (CI= red, GI*= green, CI** = 
blue). It is apparent from this figure that 1) in some cases 
the C I fine structure transitions fall so close to each other 
that they are blended, and 2) the multi-component ve- 
locity structure found in most DLAs further complicates 
the analysis. As denoted in Figure [IJ three C I velocity 
components are required to fit this DLA. Component 1 
at z a f, s —2. 625808, or v ~ —43 km s _1 , with Doppler 
parameter of b = 3.25km s _1 , is the weakest compo- 
nent and contains the largest errors. Component 2 at 
z a b s — 2. 6263247 with b = 2.57km s _1 , is located at v — 
km s _1 , and has an upper limit on N(C I**). Component 
3 at v ~ +14 km s _1 , or z a k K =2. 626491, is t he narrow 
component reported in IJorgenson et al.1 (|2009l ) with b = 
0.33 km s- 1 . 

While the velocity structure of CI in DLA 0812+32 is 
well constrained by the fitting of many multiplets, our 
confidence in the fit is increased by the similar velocity 
structure seen in neutral chlorine, CII A1347. Because 
of their similar ionization potentials, both below a Ryd- 
berg, CI and CII A1347 are often obs erved to have the 
same velocity structure pural Il974al ) . This is seen in 
Figure [2] along with several other sub-Rydberg and low- 
ion transitions that trace the velocity structure of C I, 
namely, Ge II A1237, Mg I A1827, Si I A 1845 and Zn II 

8 The abundance ratio with respect to solar is defined as [X/Y] 
= log(X/Y) - log(X/Y) 



A 2062. 

This DLA also contains relatively strong H2, with a 
total log N(H2) = 19.90 cm~ 2 , giving a molecular fraction 

of f = 0.067, where f is defined as f = n(hi)+2N(h 2 ) • ^ e 
H2 velocity components are consistent with those of C I. 

•DLA 0812 + 32, z abs = 2.066780 : Data for this 
DLA is the same as that for DLA 0812+32, z aos — 
2.62633 discussed above. There is one CI velocity com- 
ponent in this DLA. To fit CI, we used only the A1656 
and A1560 multiplets because of heavy blending with the 
Lyman a forest blueward of restframe ~A 1467 A, see 
Figure [31 This DLA also contains CII, however it is in- 
teresting to note that the centroid of the CII profile is 
displaced by «+5km s _1 with respect to v = km s _1 
at the C I centroid of z a b s = 2.066780. While this is unex- 
pected, it is also true of all of the low- ions, see Figure |H 
However, Si I and Mg I appear to align more closely with 
CI. 

•DLA1331 + 17 : IWolfe & Da'visl (|1979ft discovered 
21 cm absorption towards QS01331+17at z aos — 1.77642, 
a similar redshift at which a DL A had previously bee n 
discovered at optical wavelengths (jCarswell et al.lll975l ) . 
They deduced a spin temperature of T s « 770 — 980 K 
by combining the 21 cm line equivalent width with the 
N(H I) obtained from Lyman a absorption. 

One of the first attempts to measure the CMB 
te mperature at high redshift using C I was made 
by IMever et al.l (fl98l using DLA 1331+17. They 
measured CI and put an upper limit on the ratio 
n(CI*)/n(CI) that allowed them put an upper limit on 
the temperat ure of the CMB , Tqm b <16 K at z a bs — 
1.776. Later. ISongalia et al. I (|1994l ) used the Keck tele- 
scope to make more precise measurements. They suc- 
ceeded in measuring the CI* transition and derived a 
CMB temperature of T = 7.4 K that a greed well w i th the 
theoretical prediction. More recently. ICui et al~l (|2005D 
discovered an unusually high level of molecular hydro- 
gen (H 2 ) in DLA 1331+17 using Hubble Space Telescope. 
They detect a molecular fraction of 5.6%±0.7% in a com- 
ponent at z a bs — 1. 776553. They derive an excitation tem- 
perature of the H2- bearing com ponent of T ex = 152 K. 
Recently ICarswell et alj (|2010() discovered a narrow ve- 
locity component of CI at z a0 s = 1.776525, that requires 

gas With Tthermai < 218 K (la). 

To maximize the UV coverage for DLA 1331+17, we 
used a Keck HIRES spectrum with instrumental reso- 
lution FWHM = 6.25 km s _1 in combination with a 
bluer UVES spectrum of resolution FWHM = 7.0 km s _1 
kindly provided by R. F. Carswell. This provided cov- 
erage of the maximum number of C I multiplets, down 
to the restframe CI A1277 multiplet at observed wave- 
length A b s « 3546 A. We included in the fit the following 
multiplets: A1656, A1560, A1328, A1280, and A1277. 

As shown in Figure[5j DLA 1331+17 requires three CI 
velocity components. Component 1, at v ^0 km s _1 , is 
the strongest component in terms of column density (see 
Table [2]) with b = 5.08 ± 0.24 km s" 1 . Component 
2, at v ^17 km s" 1 , has a narrow velocity structure (b 
= 0.55 ± 0.13 kms" 1 ), and is discussed in detail in 
ICarswell et alJ (|2010l ). The third component, at v ~20 
km s _1 , with b = 24.45 ± 6.2 km s _1 , does not ex- 
hibit measurable CI* or CI** absorption. Lending con- 
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fidence to the fit is the presence of CI I, shown in Figure|5] 
along with several other sub-Rydberg and low-ion tran- 
sitions that trace the velocity structure of C I, namely, P 
II A1152, S II A1250, Mn II A2576, and Mg I A2852. 

We can compare our results to another re- 
cent measurement of CI absor ption in this system. 
iDessauges-Zavadskv et al. I (|2004f ) used VLT and Keck 
data to measure CI in two components, log N(CI) = 
13.12 ± 0.02 cm" 2 at z abs = 1.776365 and log N(CI) = 
12.72 ± 0.02 cm" 2 at z abs = 1.776523. While in the lat- 
ter case, what we call component 2, their measurements 
agree with ours to within la errors, the component 1 is 
in disagreement, with small errors, by 0.12 dex. How- 
ever this difference is reasonable considering continuum 
placement uncertainties and the fact that they did not 
include fitting of the C I fine structure states which in- 
volve considerable blending with the resonance line. 

•DLA 1755 + 578 : DLA 1755+578, at z abs = 1.9692, 
contains 8 C I components, making it one of the more 
complex C I systems. It is interesting in its own right be- 
cause we have discovered the presence of Si II* absorption 
in this DLA. While Sill* has been observed in GRB- 
DLAs pu mped by the UV rad i ation field of the GRB 
afterglow (jSavagho fc Falll 120041: IProchaska et al.ll2006h . 
this is the first known case of Sill* observed in a high 
redshift QSO-DLA. Because Sill* is generally thought to 
arise in warm gas - its excitation energy is 413K - this 
system is not only unique but also extremely interesting 
since it appears that several velocity components contain 
both C I and Si II* absorption, indicating that this is one 
of the rare DLAs exhibiting both CNM and WNM gas 
at high- z. We discuss this object in greater detail in a 
future paper. 

•DLA 2100 - 06 : While this DLA, at z abs = 3.09237, 
exhibits measurable C I, our data are not of sufficiently 
high quality to obtain good measurements of the C I fine 
structure states. This is in part because the CI lines 
are not strong, and in part because the strongest multi- 
plet, A1656, falls close to an order gap and as a result is 
contained in a lower signal to noise region. 

All CI multiplets A1277 and redwards are included in 
the fit. The fit requires three CI velocity components, 
only one of which has measurable N(CI**), albeit with 
large enough errors that it might best be considered as 
an upper limit. We present spectral coverage of the mul- 
tiplets used in Figure [7] 

•DLA 223 1 - 00 : DLA 2231- at z abs = 2.066 was 
analyzed by iProchaska fc Wolfel (11999D. who rep orted 
the detection of titanium, and " IProchaska 1 (1999) ana- 
lyzed a Lyman-limit system at z a f, s =2.652 in this sight- 
line. With an absorption redshift of z abs = 2.066 in the 
line of sight to a quasar with emission redshift of z em = 
3.02, both the A1328 and the A1560 C I multiplets of this 
DLA fall in the Lyman-a forest. While most of the A1560 
multiplet is unblended, the A1328 multiplet contains a 
blend with a large Lyman-alpha forest line. Our analysis 
includes only the A1656 and A1560 multiplets, which are 
shown in Figure |U The best fit required two C I velocity 
components: component 1 at v = —77 km s _1 , or z abs 
= 2.06534, with no measurable fine structure lines, and 
component 2 at v = km s _1 , or z abs — 2.066122, with 
detected C I fine structure lines. These two C I compo- 
nents trace the two strongest components of the other 



low ions, as seen in Figure While CI I is located in a 
region of low S/N, the detection of Mg I A 2026 lends 
confidence to the CI detection. We note that while the 
large error of the N(CI**) measurement of component 2 
indicates that it should probably be considered an up- 
per limit, for the purposes of the present analysis we will 
treat the value as a detection. 

•DLA 2340 - 00, z abs = 2.054924 : The complex 
multi-component nature of the C I in this z abs = 2.054924 
DLA made it difficult to fit despite the fact that seven C I 
multiplets (A1270, A1276, A1277, A1279, A1328, A1560 
and A1656) have been included in the fit. The A1279 
multiplet was only partially used due to blending with 
interlopers. Nine C I components were required in the 
fit, see Figure [TO] The components are labeled 1 — 9 
and located at the following velocities relative to an arbi- 
trarily chosen z abs = 2.054151 for v = km s _1 : v ~ 13 
km s _1 , 37 km s _1 , 44 km s _1 , 52 km s _1 , 55 km s _1 
, 57 km s _1 , 83 km s" 1 , and 96 km s _1 . 

DLA 2340—00 also contains a relatively high total col- 
umn density of molecular hydrogen, logN(H2) = 18.20 
cm~ 2 , giving an H2 fraction f — 0.014, making it one of 
the few relatively H2- rich DLAs. 

•ABSL 2340 - 00,z abs = 1.36 : The sightline to- 
wards J2340— 00 contains a second set of CI absorption 
lines that require 4 CI velocity components to fit. Only 
the A1656 and A1560 CI multiplets were available to 
constrain the fit because of the low redshift. Because 
of blending, in the case of A1656 multiplet with the C I 
A1279/A1280 multiplet of the high- z absorber, and in 
the case of A1560 multiplet with the Lyman a forest, the 
C I lines were fit by tying them to their associated Mg I 
transitions. The results are given in Table [2] 

3. APPLYING THE CI FINE STRUCTURE TECHNIQUE 

In this section we detail the processes by which we 
determined the constraints on the physical conditions in 
each DLA. 

3.1. The Steady State Equation 

We developed an in- house code based on POPRATIO 
(jSilva &: Viegas 1 120011 ) to calculate the theoretical C I 
fine structure level populations by maki ng the standard 
assum ption of steady state. Following iSilva fc Viegas I 
(|2001f ) we considered spontaneous radiative decay, direct 
excitation by the CMB, UV pumping, and collisional ex- 
citation and de-excitation. As in POPRATIO, the rate 
equations leading to steady state populations of state i 
is given by 

Ej n 3 i A n + b 3* u j* + r> + Ek nk <iji ) = m 

n t Ei {Aij + B ijUij + r tf + J2 k n k q% ) 1 > 

where Ay are the spontaneous decay transition proba- 
bilities, By are the Einstein coefficients, u, 3 is the spec- 
tral energy density of the radiation field, Tij is the indi- 
rect excitation rate du e to fluorescence and is defined by 
ISilva fc Viegas I (|2001| ). The quantity n k is the volume 
density of the collision partner k, where k — (H°, n e , 
n p ) and q*j- = < av >, is the collision rate coefficient. 
The reverse rates are calculated using the assumption of 
detailed balance. All coefficients were taken to be the 
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same as those used in POPRATIO, with the exception 
of collisions with neutral hydrogen, for which we used 
the more recent ra te coefficients calculated usi ng the an- 
alytical formula by lAbrahamsson et al. I (2007) extended 
to the temperature range of T = 10,000 K. For the sake 
of simplicity, we did not consider excitation by collisions 
with either molecular hydrogen (H2) or helium (He I). 
In the case of the former, the paucity of H2 found in 
DLAs, at typical fractions of less than s» 10~ 5 renders 
the effect of collisions with H2 so small as to be negli- 
gible. However, even in DLAs in which the H2 fraction 
is relatively large (for DLA 1331+17, the molecular frac- 
tion was determined bv lCui et al. I (|2005l ) to be f = 0.056 
or 5.6% ± 0.7%), H 2 does not have a large effect on the 
collisional excitation of C I. 9 In the latter case of He I, 
collision rates are significantly lower than those of other 
partners. Additionally, the density of He I compared 
with that of H I is typica lly n(He I) = 0.0975 n(H I) 
([Anders fc Grevesse II1989T ). making collisions with He I 
relatively unimportant. 

Direct excitation by the CMB is calculated assum- 
ing the standard cosmology and a CMB temperature o f 
T = T (l + z) where T Q = 2.725 K ([Mather et al. II1999D . 
At high redshift, the CMB radiation generally has the 
strongest effect on the C I fine structure level populations 
because of the small temperature difference between the 
ground and first fine structure states. However, depend- 
ing on the physical circumstances, other mechanisms 
such as UV pumping or collisions can dominate. 

We included UV pumping due to a radiation field con- 
sisting of two components that we will call J v Bkd , and 
J u local , and let J u total = jJ ocal + J v Bkd . J v Bkd is the 
background radiation due to the integrated contribution 
from high z galaxies and quasars, known as t he Haardt- 
Madau background (|Haardt fe Madaul (jT996) ; and more 
recently using CUBA 10 ). In all cases, the minimum value 
of the total radiation field is determined by J v Bkd . In 
each case the value of J v Bkd is calculated based on the 
redshift of the DLA, and these values are summarized 
in Table JH For an explanat ion of these values, see Fig- 
ure 1 of I Wolfe et al. I ((2004). Because each CI -bearing 
DLA also contains strong CII* A1335.7 a bsorption, we 
used t he C II* technique (see WPG03 and I Wolfe et all 
(2004)) to estimate the local radiation field due to star 
formation, j v local >cn ^ anQ i mcmc i ec i this contribution 
in the UV pumping. The jJocaLCir is ca ,l cu l a tcd at 
A = 1500A, or 8.27 eV, and in § 15X21 we explain the 
estimated error on this value, ~±50% . Note that in 
general, J v local > on * _ jjocal^ ry e introduce the nota- 
tion J u local > G 11 to specify how the local radiation field 
is measured, i.e. in this case it is the local star formation 
rate per unit area measured via the C II* technique. 

To quantitatively determine the effects of the radiation 
field on the CI fine structure level populations, we plot 

9 I.e. the rate of excitation due to collisions with neutral hy- 
drogen at T = 100 K is n H 1 q^ 1 = 2.976 xlO" 10 n(H I) s~\ 
where we use the conventional notation of representing CI, CI*, 
CI** by the indices 0, 1, 2 respectively. If we take the fraction of 
H2 to be 5.6% of H I, the rate for collisions with H2 is n H2 q^ 2 
= 6.7 xlO" 11 (0.056) n(H I) = 3.75 xl0~ 12 n(H I) s" 1 , which is 
roughly two orders of magnitude smaller than qj^ 1 . 

10 CUBA (Haardt fe Madau. 2003) is available at: 
http: / / pitto.mib.infn.it / ~haardt / cosmology.html 
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for excitation by a wide range of radiation fields in Fig- 
ure 111! In this example we have considered collisions 
with neutral hydrogen at T = 100 K, spontaneous ra- 
diative decay, and the excitation by the CMB at z = 2, 
in addition to a radiation field of varying strengths as 
denoted in Figure [IT] At low density, excitation by the 
CMB is dominant, unless the input UV radiation field 
is strong, in which case, the UV dominates. At higher 
densities, i.e. n(H I) > 10 cm -3 , collisional excitation be- 
comes important and finally at n(H I) >1000 cm~ 3 the 
levels thermalize. In other words, at low densities the 



CMB sets the floor of 



n{C I) ■ 



Only when the radiation 



field exceeds a total strength of approximately J, 



total 



> 10" 18 - 5 ergs cm" 2 s" 1 Hz" 1 sr" 1 , 



does be S in 



to significantly exceed that caused by the CMB alone. 
Since J v Bkd rarely exceeds 10 -19 ergs cm~ 2 s _1 Hz -1 
sr _1 , the Haardt-Madau background alone at a redshift 
of z ~ 2, has essentially no effect on the C I fine structure 
excitation. 

3.2. Steady State Solution 

The steady state densities in each of the C I fine struc- 
ture states are found by solving three homogeneous equa- 
tions with three unknowns. To solve the three homoge- 
neous equations, we solve for the ratio of each excited 
state relative to the ground state. We denote all terms 
involving CMB excitation, UV pumping and collisions by 
the shorthand, Rjj = B^- u^ + + ^2 k n qfj which is 
summed over the k different collision partners, where all 
terms were defined in the previous subsection. Reverse 
reaction rates are calculate d through the prin c iple o f de- 
tailed balance. Following Jen kins fc Tripp 1 (|200lD (see 
their equations 10 and 11), we find: 



n(cr) _ 

n(CI) - (9] 

(Roi)(A 2 , 1+A 2 ,0+B. 2 . l+fi2 0) + (fld,2)(A21+fl2,l) \ A > 

{Ri,2){A 2sl +R2,a) + {A 1 ,o+Ri,tt)(A 2 ,i+A 2sl +R2,i+R2,o) 

and 



n(C7") _ 
n(C I) ~ 

(flo,2)(A 2 ,0+-R.l,0+-Bl,2) + (.Ro,l)(-Rl,2) 

{Ri,2)(A 2:0 +R2,a) + (A 1 ,a+R lfi )(A2,i+A 2 , +R2,i+R2,a) 



(3) 



where the states C I, C I*, C I** are denoted by the indices 
0, 1, 2 respectively. The resulting theoretical solutions 
are functions of the density of neutr al hydrogen, n(H 
I), an d the temperature. Following Uenkins fc Shaval 
(|1979D we define, 



n(CJ* 



n(C I*) 
n(C I) 



and 



/2 



n(CI) t 



n(CI* 
n(CI)t 



1.0 



1.0 



n(CI*) , n(CI") 
n(C I) + n(CI) 



n(CI") 

n(C I) 

n(C I*) , n(C /") 
n(C I) 



(4) 



(5) 



MCI) 



wheren(C I) to t = n(C I) + n(C /*) + n(C /** ). We give 
the values of (fl, f2) for each component of each DLA in 
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Table p 

In Figure[12]we plot the theoretical solutions in the (fl, 
f2) plane for the example case of component 3 of DLA 
0812+32. In this case, jJ°tai,Cir as derived from the 
CII* technique, is j^taLCir = 7 4 xl0 -i9 er gs cm" 2 
s _1 Hz -1 sr -1 . The solutions follow tracks, one for each 
temperature from T = 10 K — 10,000 K increasing in 
steps of 0.1 dex. Along each track, n(H I) increases from 
10~ 3 - 5 cm" 3 to 10 41 cm" 3 in steps of 0.02 dex, each den- 
sity being represented by an individual point. We chose 
the ranges of temperatures and densities to ensure that 
we had broad coverage of the entire plane of possible so- 
lutions. While the C I data may allow higher T values, 
these would imply a primarily ionized gas which is highly 
unlikely for this material. The data point at (fl, £2) sa 
(0.31, 0.09) is determined by our curve of growth anal- 
ysis (in the case of the resonance line CI) and VP FIT 
fits to the data for component 3 (see Tabled]), la (2a) 
errors are determined by calculating the values of fl and 
f2 using the upper and lower allowed values as deter- 
mined by the la (2a) error bars for each column density 
measurement. The resulting red (blue) polygon defines 
a region of space in the fl versus f2 plane that contains 
the acceptable la (2a) solutions. In Figure Q2] we plot 
the predicted theoretical solutions that fall within the 
la (in red) and 2a (in blue) error polygons of the data, 
on a graph of density versus temperature. These are the 
allowed density and temperature combinations derived 
from the C I data under the assumption of a single phase 
absorber. 

3.3. Ionization Equilibrium 

To further constrain the C I solutions we introduce the 
assumption of ionization equilibrium and utilize the mea- 
sured ^tj- ratio, which we assume is equal to = 

^j^jjy. Ionization equilibrium is a reasonable assump- 
tion because the densities in these clouds are generally 
relatively high, thereby ensuring that the recombination 
times are shorter than the typical dynamical timescales. 
In basic form, ionization equilibrium can be written: n e 
n(X + ) a = n(X) T, where a is the case A recombination 
coefficient of element X + to X and T represents the ion- 
ization rate. From this equation we get, ^~rgrf^ = ^-j, 

i.e. the ratio ^jj- is a function of T (where T is propor- 
tional to the radiation intensity) , n e , the electron density, 
and a, which is a function of temperature. Therefore, for 
a given radiation field and temperature, we determine 
^jj- for a range of possible electron densities, n e (De- 
tails of the ionization equilibrium are given in WPG03, 
and we discuss the sensitivity to the radiation field in 
section 13.3. 1|) . We then use the measured ^tj- of each 
DLA to constrain our allowed solutions. While we mea- 
sure N(CI) directly, we must employ the conventional 
method of measuring N(C II) by proxy using N(Si II), be- 
cause the available CII A1334 transition is saturated in 
all the DLAs considered here. We generally measure the 
Si II or S II and Fe II using the standard AODM, which is 
well suited for cases such as these in which we have more 
than one transition of an ion. As in I Wolfe et all (2004) 
we let [C/H] = [Si/H] + [Fe/Si] mt where the intrinsic 
(nucleosynthetic) ratio [Fe/Si]i„t = —0.2 for a minimal 



depletion model or [Fe/Sil mt = 0.0 for a maxima l de- 
pletion model. We follow iMurphv fc Liskd (|2004ft and 
adopt the minimal depletion model in this work, and an- 
alyze the implications of the minima l versus maximal de- 
pletion model in section l5XTl While I Wolfe et al. I ()2004f ) 
estimate the error in the measurement of ^ff to be 0.1 
dex, we use the more conservative estimate of 0.2 dex. 
We also note that, assuming there is no hidden satura- 
tion of metals, this is the maximum N(CII) that could 
be associated with the C I gas, and we discuss the impli- 
cations of this assumption further in section 15.1.31 

To demonstrate the constraints imposed by ionization 
equilibrium, in Figure [T3] we plot the C I solutions in 
terms of n(H I) versus log(^y-) for the example case of 
component 3 in DLA 0812+32. The measured log(%y-) 
= 3.10 is indicated by the red dashed line with a range of 
±0.2 dex indicated by green dashed lines. It is obvious 
that the region of allowed CI fine structure solutions is 
further constrained by invoking ionization equilibrium. 
For clarification, in Figure [l_5]we re-plot the density ver- 
sus temperature diagram, however now we denote the 
final solutions, those constrained by the ^rj- ratio, in 
black (lcr) and yellow (2a). As a result, invoking ion- 
ization equilibrium results in even tighter constraints on 
the densities and temperatures of the C I -bearing cloud 
without making any assumptions about the gas temper- 
ature. In this example case we constrain the density to 
be 72 < n(H I) < 549 cm" 3 while the temperature is 
constrained to be 25 K < T < 251 K, and the pressure 
3.54 < log (P/k) < 4.80 cm -3 K. x 2 minimization finds 
the best-fitting solution: n(H I) = 100 cm" 3 , T = 79 K, 
and log (P/k) = 3.90 cm" 3 K. 

3.3.1. Using Ionization Equilibrium to Constrain the 
Radiation Field 

Until now we have assumed the UV radiation field due 
to star formation, J v local , is determined by the CII* 
technique. We will now relax this constraint and al- 
low the input local radiation field to vary, repeat the 
above analysis for each case, and examine the results as 
a function of radiation field. We show that for some C I - 
bearing clouds we can place upper and lower limits on the 
allowed radiation field utilizing only the C I fine structure 
data and the assumption of ionization equilibrium. 

We use a grid of Jj ocal values from J v local w 10" 21 - 
10~ 16 ergs cm" 2 s" 1 Hz" 1 sr" 1 . To differentiate these 
radiation fields from those predicted by the C II* tech- 
nique, we will call them j v local fi 1 . For each j v local > c 1 wc 
first add the J u Bkd to obtain a J u totaLC 1 , and then re- 
run the above analysis, calculating the fl versus f2 tracks 
for each j v total c 1 ; followed by the ionization equilibrium 
analysis. 

In Figure [TBI we plot the lcr CI results for the example 
case of DLA 1331+17 on a graph of ^Jj- versus n(H I) 
where for clarity we show only a sub-sample of the entire 
jjotai,c i g r j^ resu it s . Solutions corresponding to each 
jjotai,c i are re p reS ented by different colors. Again, we 
use the ^J-j- data to constrain the allowed solutions, and 
it is apparent that for j v total & 1 > \ xl0" 19 ergs cm" 2 
s" 1 Hz" 1 sr" 1 , there are no acceptable solutions. This 
limit on j u total fi ' 1 j s actually quite strict considering that 
the Haardt-Madau background is only ~1 order of mag- 
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nitude lower than this value. 

This technique places an upper limit on the local ra- 
diation field, and hence the star formation rate, using 
only C I fine structure absorption and the assumption 
of ionization equilibrium, without invoking any assump- 
tions about star formation from the C II* technique. We 
note in this example, of DLA 1331+17, component 3, 
that the strength of the radiation field as determined by 
the C II* technique (JJotal,Cii* _ 3 3 xl0 -i9 er gs cm" 2 
s _1 Hz -1 sr" 1 ) is just slightly higher than the la upper 
limit derived independently from the C I data. However, 
in the case of DLA 1331+17, the C II* transition is likely 
blended with a Lyman— a forest line, forcing us to esti- 
mate the true N(C II*) by assuming it traces the velocity 
structure of other low ions and fitting it together with 
the forest line as explained in section 14.31 

A byproduct of this analysis places limits on the den- 
sity, temperature and pressure of the cloud. In this ex- 
ample case of DLA 1331+17, at the la level the density 
is limited to the range, 11 < n(H I) < 44 cm" 3 while the 
temperature is constrained to be 79 < T < 794 K, and 
the pressure 3.50 < log (P/k) < 4.04 cm -3 K. Again, 
we stress that these limits are derived independently of 
the results of the C II* technique, using only the C I fine 
structure data and the assumption of ionization equilib- 
rium for a range of possible radiation fields. 

4. RESULTS 

A summary of the details of each DLA are presented 
in the following tables: Table [2] a summary of the C I 
data, Table [3J a summary of general DLA information, 
Table 21 the resultant radiation fields derived from the 
C II* technique, and in Table [5] the final C I fine struc- 
ture solutions giving allowed ranges of densities, temper- 
atures and pressures. Finally, in Table [Bl we give the 
results of lifting the C II* constraint on the local radia- 
tion field. 

We now briefly describe the results for each C I -bearing 
DLA. The casual reader may wish to skip directly to the 
more general discussion of the results in § [5j 

4.1. FJ0812+32, z abs =2.62633 

This DLA contains three distinct CI velocity compo- 
nents whose CI fine structure levels we will analyse in- 
dividually. In measuring the metallicity, N(C II) and 
N(CII*) we will take two approaches; The first, dis- 
cussed here, and the second, the individual component 
analysis, discussed in section 14.1.11 Heavy blending of 
line profiles and the uncertainty / degeneracy of line pro- 
file fitting techniques have led to the customary use of 
the AODM technique to measure the amount of metals 
over the entire, blended DLA profile. Additionally, this 
avoids any question about the distribution of the N(H 
I), which by definition is damped and therefore kinemat- 
ically unknowable. Because we are measuring the metals 
over the entire (blended) DLA, we refer to these mea- 
surements as 'global', and we apply them to each CI 
component, which we consider reasonable given that the 
low-ions such as S II and Si II track the velocity compo- 
nents of C I, see Figure [5] 

We measure log N(CLT) = 14.30+0.01 cm' 2 which 
along with the C II* technique determines that the radi- 
ation due to stars jJ° cal > CH * = 7.17x 10" 19 ergs cm" 2 



s 1 Hz 1 sr . Adding in the contribution from the 
background, J v Bkd = 2.45 x 10" 20 ergs cm" 2 s" 1 Hz -1 
sr -1 , gives Jjota^cw ^7.4 x 10" 19 ergs cm" 2 s" 1 Hz -1 
sr -1 . The excitation rates due to the local star formation 
alone are r i = 3.17x lO^s" 1 , r 02 = 2.40 XlO'V 1 , 
and Ti2 = 3.12 xl0 _9 s _1 . These can be compared with 
those of the Haardt-Madau background at the redshift 
of this DLA that produces excitation rates Tqi = 1.09 

xio- 10 s-\ f 02 = s^xio- 1 ^- 1 , r 12 = 1.07xl0- 10 s- 1 . 

We used the total radiation field and the assumption 
of ionization equilibrium to determine solutions for the 
three CI components that are labeled 0812+32 g i o b a z in 
Table [5l We are unable to perform an analysis of com- 
ponent 2 given that the upper limit of N(CI**) < 12.39 
cm" 2 is relatively large, resulting in a (fl, f2) combina- 
tion with no limitation on n and T. 

We can use the constraints on the volume density to 
estimate the physical size of the cloud. Using the neutral 
hydrogen column density logN(H I) = 21.35 cm" 2 and, 
for example, the best-fit volume density for component 
3, n(H I)« 100 cm" 3 , we can estimate the size of the CI 
-bearing cloud to be I = N(H I)/n(H I) k. 2.24 xlO 19 
cm, or « 7 pc. This estimation is technically an upper 
limit to the size of the cloud that assumes that all of the 
H I is associated with the CI component. 

When we relax the constraint of J v total > c 11 as derived 
from N(C II*), as discussed in the previous section, at the 
2 a level we constrain the density to n > 0.1 cm" 3 and n 
> 7 cm -3 for components 1 and 3 respectively. In both 
cases we place a not-so-strict upper limit on the allowed 
radiation field of J u total < 773 xlO" 19 ergs cm" 2 s" 1 
Hz" 1 sr" 1 . 

4.1.1. DLA 0812+32: Individual Component Analysis 

While the initial modeling was completed using the ra- 
diation field and metallicity as measured over the entire 
profile, it is clear upon careful i nspection of the spectru m 
(see Figure [2] and Figure 7 of Uorgenson et al . (2009)), 
that the depletion, and hence, dust-to-gas ratio, is not 
constant over the three components. This motivates our 
attempt to analyze each velocity component individu- 
ally. Specifically, component 3, the narrow component, 
contains obvious Zn II, C II* and C I, while it exhibits no 
evidence for Cr II. A natural explanation is that Cr II is 
heavily depleted onto dust grains in compone nt 3. How- 
ever, as discussed in Uorgenson et all (|2009f ). blending 
precludes a straightforward measurement of the equiva- 
lent width, and profile fitting using VPFIT, produces a 
Zn II column density in component 3 that is unrealis- 
tically high and ruled out by an upper limit on N(0 I) 
assuming solar relative abundances. We explain this by 
assuming the prese nce of another weak, b road compo- 
nent (for details see Uorgenson et al.l (|2009D ). and model 
the log N(Zn II) in component 3. 

Because Zn II and Cr II are not saturated and ap- 
pear to trace the velocity structure of C I quite well, we 
use these ions to determine the metal distribution, the 
dust-to-gas ratio, and the N(C II) in each component 
separately. To do this we first tie the redshifts and b 
values of each ion together and then use VPFIT to de- 
termine the column densities in each component. It is 
expected that Zn II and Cr II be tied together because 
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they are metal-line transitions arising in singly ionized 
species with similar ionization potentials, and thus are 
likely to show the same physical/ velocity structure. How- 
ever, they will not necessarily exactly trace the C I, given 
that CI is affected by the incident J u total and n e , and if 
there is a gradient in J l/ total or n e , this could cause a 
difference between the structure of C I and the other low 
ions. We give the re sults of this fitting in Table [7] In 
light of the work of Uorgenson et al.l (|2009() we include 
an additional broad weak component, in order to achieve 
a realistic N(Zn II) in component 3. 

Due to the nature of damped lines, we cannot use 
Lyman-a to determine the distribution of the neutral 
hydrogen among these components. Therefore, we as- 
sume that the N(H I) traces the low-ion metals, in this 
case Zn II, and that the neutral hydrogen is distributed 
proportionally to the metals. This results in each com- 
ponent having the same metallicity, [Zn/H] = —0.58. We 
note, that in theory, this metallicity should be the same 
as that determined for the global case, [Zn/H] = —0.81. 
However, the individual component analysis reveals the 
presence of the narrow component 3, for which we have 
estimated the N(Zn II) based upon an upper limit on 
N(0 I). While the exact metallicity of this component 
remains a mystery without higher resolution data, it is 
apparent that component 3 lacks significant Cr II, in- 
dicating a high level of depletion, and consequently a 
higher dust-to-gas ratio than components 1 and 2. Fol- 
lowing WPG03 (their equation 7) we define the dust-to- 
gas ratio relative to the Milky Way as follows, 



K _ 1Q [Zn/H\ int f 10 [Cr/Zn] int _ 1Q [Cr/Zn] ga3 \ ^ 

where "gas" is the abundance ratio in the gas phase and 
"int" is the intrinsic nucleosynthetic abundance. We cal- 
culate the dust-to-gas ratio in component 3 to be ^17% 
of the Milky Way (Ior k = —0.78). For comparison, 
Prochaska_et_aLJ (j2003D measure k ~ 6% (logft=— 1.24) 
over the entire profile of this DLA, and the typical dust- 
to-gas ratio in DLAs is ^1/30 solar. We summarize in- 
formation about each component in Table [7] 

As we might expect, given that the column densities 
of all low ions with the exception of Cr II, but includ- 
ing C II* , are much higher in component 3 (due to the 
fixed small Doppler parameter), the j v local > c 11 of com- 
ponent 3 deduced from the C II* technique is more than 
1 order of magnitude larger than that determined for 
components 1 and 2. The results of the CI models for 
each component are presented in Table [5] where they are 
denoted by the subscript 'div'. Note that the radiation 
field for component 3 predicted by the CII* technique 
is quite large, at 260 x 10~ 19 ergs cm~ 2 s _1 Hz -1 sr _1 . 
The solutions require low temperatures (T < 32 K) and 
densities in excess of 10 3 cm -3 , which imply an upper 
limit on the cloud size of I < 0.1 pc. This small size may 
be in conflict with the evidence against partial covering 
provided by the H 2 , see section 14.1.21 It is possible that 
we have over-estimated the amount of N(CII*) associ- 
ated with the narrow component 3, and that more of the 
N(C II*) is associated with the broader component 4, but 
we would require a higher resolution spectrum to confirm 
this scenario. 

Relaxing the constraint of J u local determined by the 



C II* technique, gives no limits on temperature, and a 
large range of allowed densities for component 1 ( n(H 
I) = 0.002 - 4166 cm" 3 ), with a lower limit of n(H I) 
> 6 cm -3 for component 3. Constraints on the allowed 
radiation fields are consistent with that determined by 
the CII* technique. 

4.1.2. DLA 0812+32 Molecular Hydrogen 

This DLA also shows evidence of relatively strong 
molecular hydrogen (H 2), see Figure [17] As shown in 
Uorgenson et al.l ([20090 . there is no evidence for partial 
coverage, as the H2 lines are black at line center. We 
used VPFIT to fit the H2, which required three com- 
ponents whose redshifts roughly agree with those of the 
three C I components. While the Doppler parameters do 
not agree with those of the C I- they are roughly similar 
to or smaller, whereas we would expect them to be larger 
by a factor of (6) 1 - we will refer to these H2 components 
as components 1,2, and 3, and assume that they are co- 
spatial with the C I components. This is not the first 
report of such differences i n the b- values between H2 a nd 
CI components in DLAs (jNoterdaeme et al. Ifeoofallbh . 
Because of the close match in redshift space and the 
likelihood that CI and H 2 are co-spatial, we fixed the 
well-determined Doppler parameter of component 3 to 
that of the C I (b CI = 0.33 k m s" 1 ), scaled by the re la- 
tive atomic weights as done in Uorgenson et al.1 ( 2009n 11 . 
We list the parameters of the H 2 fits in Table [H 

It is interesting to note that like the C I, the majority 
of the H2 resides in component 3, the narrow velocity 
component. The fraction of H2 in component 3, f > 
0.06, or > 6%, with the upper limit representing the un- 
certainty in N(H I) distribution, is among the highest 
found in DLAs. While we have tied the Doppler param- 
eter to that of the C I, we note that because of heavy 
saturation of the J = and J = 1 lines in this compo- 
nent, log N is determined by the damping wings of the 
profile and is therefore insensitive to the choice of b. We 
have verified this by artificially fixing the Doppler pa- 
rameter to a range of values from b = 0.2 - 6 km s _1 as 
well as by raising the continuum by 5%. None of these 
tests change the resultant N(H2, J=0) and N(H2, J=l) 
by more than ^0.03 dex. While the J = and J = 1 lines 
are completely insensitive to the Doppler parameter, the 
higher level J states are sensitive to the choice of Doppler 
parameter. Therefore, we also report a model in which 
we have allowed the Doppler parameter of component 3 
to be determined by the J = 3 state and we report the 
results, called model 2, in Table [8] 

We can estimate the kinetic temperature of the clouds 
using the column densities of H 2 in the J=0 and J=l ro- 
tational states and assuming the states are thermalized 
acc ording to the Bolt z mann distribution (see equation 8 
in iLevshakov et al.1 ((2002)). The excitation tempera- 
ture, T ex , is defined by, 



N(J) _ g(J) 
N(0) .9(0) 



exp[ 



B v J{J+l) 



(7) 



where B„ = 85.36 K for the vibrational ground state and 
g( J) is the degeneracy of level J. In Figure [IS] we show 

11 b H * = ( 6 )i/2 b Ci = q gi km s -i , 
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the standard H 2 excitation diagram, log(N/g) versus en- 
ergy, for the J=0 to J=5 levels for each component. As 
explained in Appendix 1 (on molecular hydrogen), the 
kinetic temperature of the gas is proportional to the neg- 
ative inverse of the slope determined by the J=0 and 3=1 
levels, assuming that the levels are thermalized (a good 
assumption given the densities of these clouds), i.e. T ex 
equals the kinetic temperature. For components 1,2, and 
3 we determine the following temperatures, T 91 = 102 
K, Tf x = 73 K, T° e l = 47 K, respectively. These are con- 
sistent with the temperatures derived from the C I data, 
which leads us to believe there is good correspondence 
between the CI and H2 data and that the gas probed 
here is a CNM. Also, the T°* = 47 K derived for the 
narrow component 3 is consistent with the upper limit 
on the thermal temperature of Tthermai < 78 K required 
by the curve of growth determined Doppler parameter. 

We can also use the H2 data to determine the am- 
bient/incident radiation (UV flux) field on the cloud, 
jj,otal,Hi ( us j n g t ne j— and J=4 states (see Appendix 
1 for details). The level populations above J = 1 arc 
unlikely to be thermalized since their larger Einstein A 
coefficients implies that these states are depopulated by 
spontaneous emission more rapidly than by collisional de- 
excitation; this rules out the condition of detailed balance 
required to establish thermal equilibrium. Instead, these 
states are likely to be populated by UV pumping, which 
is wha t we shall assume. Following iHirashita h, Ferrara I 
(2005) we call this radiation held J v , as determined 
by absorption in the Lyman- Werner H2 bands. Since the 
H2 measurement should reflect the total incident radia- 
tion field, J u local + J u Bkd = J V LW = j v total ^. Results 
for the photoabsorption rate, /?, and for the shielding 
terms, both self-shielding by H2 and shielding due to 
dust, are given in Table [7] Note that /?o is derived from 
the J = 4 population, while /?i is derived from the J = 
5 population. Using the J = 4 population, we obtained 
the following radiation fields for components 1,2, and 3 
respectively: JJ°^ H * = < 3.7 xlO" 20 , 1.7 xlO" 20 , 3.6 
xlO" 20 ergs cm" 2 s" 1 Hz -1 sr" 1 . Note, that for compo- 
nent 1, N(H2, J = 4) is an upper limit, and hence, this 
is technically an upper limit on the radiation field. Also 
notice that these radiation fields are only slightly larger, 
or smaller in the case of component 2, than J v Bkd . This 
places a rather strict limit on additional radiation from 
local star formation. 

Interestingly, the radiation field as calculated from the 
H2 is rather independent of the choice of global versus 
individual component parameters. In the above we gave 
the results of the individual component analysis, where 
we assumed that the N(H I) tracks the metals, in this 
case the N(Zn II), and used the individual component 
metallicity, N(H I) and dust-to-gas ratio to calculate the 
radiation field. If we instead use the 'global' model, the 
N(H I) increases to the global value, here log N(H I) = 
21.35 cm -2 , while k decreases in each component to the 
global value, here k = 0.07. These two effects essentially 
cancel each other out in the calculation of shielding (see 
equation [501 in Appendix A) and we obtain very sim- 
ilar radiation fields in each case (i.e. for component 1, 
J v totalH2 global < 4.2 xlO" 20 ergs cm" 2 s" 1 Hz" 1 sr" 1 , 
compared with J v total ' H2 indiv < 3.7 xlO" 20 ergs cm -2 
s" 1 Hz" 1 sr" 1 ). 



We do not understand why the radiation field derived 
from the H2 data is ~ 20 times smaller than that derived 
by the CIP technique, J v total = 7.4 xlO" 19 ergs cm" 2 
s _1 Hz" 1 sr" 1 . Additionally, the radiation field predicted 
by the H2 levels for component 2 is low enough that it 
is excluded by the C I data and the Haardt-Madau back- 
ground (however, we note that the difference is not large 
and it could be within the errors - i.e. for component 2, 
if we consider the errors on the J = and J = 4 levels 
we derive J v total ^ = 3.1 xlO" 20 ergs cm" 2 s" 1 Hz" 1 
sr" 1 , consistent with the Haardt-Madau background, but 
still « 20 times smaller than the CIP prediction). This 
conflict could potentially be resolved if either 1) we are 
underestimating the population in the N(H2, J=4) level 
- which could be caused by hidden saturation of narrow 
components as was demon strated for the narro w com- 
ponent in DLA 1331+17 (jCarswell et all 120101) . where 
the N(H2, J = 4) is likely underestimated because of the 
presence of the small Doppler component which contains 
the majority of th e molecular gas (see f urther discussion 
in section Ep] and iCarswell et ail (|2010D 1 or 2) the CIP 
prediction is overestimating the radiation field. In the 
latter case, the radiation field is dependent upon the as- 
sumption of the equilibrium pressure existing at the ge- 
ometric mean (see Section 15.1.11 for details). If instead 
the equilibrium pressure is located at P m ax or P m j n , the 
resultant radiation field can change by up to a factor of 
10 (for example, see Table \V2l) . In the former case, we 
estimate that if the radiation field as predicted by the 
CIP technique is accurate and the H2 feels this entire 
radiation field, we require the amount of H2 in the J=4 
level to increase by a factor of 20, to N(H2, J=4) ~ 15.3 
cm" 2 . This large amount of H 2 is ruled out by the data 
unless it exists in a very narrow component with b < 0.2 
km s" 1 . 

We also derive densities from the H 2 data (see Ap- 
pendix 1), assuming the temperature is equal to the ki- 
netic temperature derived from the J = and J =1 states. 
The resultant densities, n(H I) = 21, 11, and 37 cm" 3 
for components 1, 2, and 3 respectively, are nearly con- 
sistent with the C I results (where for components 1 and 
3 we derived n(H I) = 23 - 417 cm" 3 and 72 - 549 cm" 3 , 
respectively). Again, since the radiation fields were not 
consistent and they are involved in this calculation, we 
may expect some of the discrepancy. We note however, 
that the densities are consistent with the lower limits de- 
rived from the CI data alone (n(H I) >0.1 and >7 cm -3 
respectively) . 

Additionally, we can use the J = 2 rotational state 
to place an upper limit on the density if the excitation 
temperature, T^ 2 is not equal to that of T 91 , indicat- 
ing that the 3 = 2 state is not thermalized, and that 
the density is below n cr u, the critical density needed for 
thermalization. We plot critical density as a function of 
temperature for the H2 J states in Figure Q1JJ In com- 
ponent 1 of DLA0812+32, T° 2 = 132 K is greater than 
T01 = 102 K, indicating that the density must be less 
than n cr i t . In this case we derive a upper limit on the 
density of component 1 of n <200 cm" 3 , consistent with 
the derived densities for this component. In component 
3, T° 2 = 57 K is greater than T i = 47 K, also indicat- 
ing an upper limit on the density of component 3 of n 
<700 cm" 3 . For component 2, T02 is slightly less than 
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Toi - indicating either that it is thermalized, or that this 
is within the errors. In either case, the upper limit de- 
rived from the next higher state, T03, would not be that 
restrictive as n cr i t ~ 10 3 cm" 3 . 

4.2. FJ0812+32, z abs =2.066780 

The C II* transition of this DLA falls in the Lya forest, 
and because the profile is quite different from that of the 
other low low ions (see Figure 01 , it is difficult to make a 
definitive estimate of the true CII* column density. In- 
stead, we attempt to estimate the star formation rate by 
assuming three mo dels motivated by the bimodality of 
DLA cooling rates (| Wolfe et al J 120081 ). In the first case, 
case (a), we assume this is a 'low cool' DLA and use the 
average low-cool l c = 10~ 274 ergs s" 1 H" 1 . This assump- 
tion is likely the closest to the truth given that several 
physical traits of this DLA match those of the low-cool 
population of DLAs, i.e. the small value of the low-ion 
velocity, A v = 26 km s" 1 is more likely to be drawn 
from the low-cool sample, with median A v — 46 km s" 1 
, than from the high cool sample with median A v = 104 
km s" 1 . The Sill 1526 rest-frame equivalent width is 
also small at, W1526 = 0.22 ± 0.01 A, compared with the 
low-cool median W1526 = 0.26 ± 0.09 A, while the metal- 
licity, [M/H] = -1.38 ± 0.01, is slightly higher than the 
median metallicity for the low-cool population [M/H] = 
— 1.74 ± 0.19. This discrepancy in metallicity may sim- 
ply be the result of metallicity increasing with decreasing 
redshift, as the z aos = 2.06 of this DL A is lower than th e 
median z abs = 2.85 of the l c sample (|Wolfe et al.ll2008l ). 
Finally, the dust-to-gas ratio logioK=— 2.74 is similar to 
the median low-cool dust-to-gas ratio, logioK=— 2.57 ± 
0.17. In the second case, case (b), we will again make 
the low cool assumption, but instead of the standard 
a-enhancement assumption of [Fe/Si] = —0.2, we will 
assume [Fc/Si] = 0.0. Finally, in the last model, case 
(c), we will assume that the DLA is a 'high cool' DLA 

and has an £ c equal to the median high cool DLA £ c = 
10 -26.6 ergs s -i H -i 

We present the results for the three cases in Table [SJ 
It is clear from the resultant densities that either cases 
(a) or (b) are more likely to be correct. In case (c), there 
are no acceptable la solutions, and the range of densities 
at the 2<T level (n(H I) > 3800 cm -3 ) seems unphysical. 
Densities this high are unlikely to be observed because 
they imply a very small cloud size that would be unlikely 
to cover the background quasar, i.e. here, the cloud size 
is estimated as t < N(H I)/n(H I) < 10 215 /3800 = 8.3 
x 10 17 cm < 0.3 pc. 

Following our assumption that the minimal depletion 
model is more likely to be correct, we conclude that case 
(a) is the best physical approximation for this system. If 
we now relax the constraint imposed by the C II* tech- 
nique, and model a range of possible j v total - cl j we find 
that at the 2a level, we can restrict the density to n > 
32 cm cm" 3 , T < 1585 K and J u total = 0.41 - 195 x 
10 _19 ergs cm" 2 s" 1 Hz -1 sr" 1 , see Table [6] 

4.3. DLA 1331+17 

Here we discuss the analysis of component 1, at z a i, s = 
1.77636, the only component with all three C I fine struc- 
ture lines detected. The Haardt-Madau background at 



the redshift of this DLA is given by J iy Bkd = 2.53x 10" 20 
ergs cm -2 s _1 Hz -1 sr -1 . To estimate the local stellar 
contribution to the UV field we utilized the CII* tech- 
nique. While WPG03 report logN(CII*) = 14.05±0.05 
cm" 2 for this object, the C II* transition is likely blended 
with a Lya forest line because its velocity profile is very 
different from that of the low ions, see Figure [6l Gen- 
erally, the C II* velocity profile traces that of low ions 
such as Sill A1808. Therefore, we cannot exclude the 
possibility that the C II* transition is blended with a for- 
est line, and we take the WPG03 value of logN(CIF) = 
14.05±0.05 cm -2 as an upper limit. To obtain an esti- 
mate of the actual N(CII*), assuming that it is blended, 
we (Carswell 2007; private communication) fixed the 
shape of the velocity profile to that of Sill A1808 and 
normalized the C II* contribution by fitting the profile si- 
multaneously with a coincident Lya line. Summing over 
the C II* components results in an estimate of the true 
N(CIF) oflogN(CII*) < 13.56 cm" 2 . With this value of 
C II* absorption, the C II* technique results in a star for- 
mation rate that produces a radiation field of J v local < c 11 
ps 3.1 x 10" 19 ergs cm" 2 s" 1 Hz" 1 sr" 1 . Therefore, the 
total input radiation field, that is, J u Bkd + jj° cal c II * 
, is J v total,cu m ~ 3 3x 1(r i9 ergs cm" 2 s" 1 Hz" 1 sr" 1 . 
Finally, the C I fine structure data constrained by ioniza- 
tion equilibrium, gives 2a results of 16 < T < 32 K, 91 
< n(H I) < 363 cm" 3 , and 3.36 cm" 3 K < log(P/k) < 
3.86 cm" 3 K . 

We can compare our results to those of iCui et al~l 
(2005), who used CLOUDY to derive a best fit model to 
their H2 data that resulted in a hydrogen number density 
n(H I) sa 0.2 cm" 3 and T ss 140 K for the H 2 bearing 
cloud at z a b s — 1.776553. This cloud is clearly not in 
pressure equilibrium with the CI cloud of our analysis, 
component 1 at z aos =l. 77637, i.e. while the tempera- 
ture is somewhat higher, the density is more than two 
orders of magnitude lower, resulting in pressures of P/k 
ps 28 cm" 3 K for the H2 -bearing cloud. The C I -bearing 
cloud more closely resembles the CI -bearin g clou ds in 
the local ISM as found by I Jenkins & Tripp I (|2007l ) and 
Uenkins fc Tripp I (|2001l ). We discuss these similarities 
f urther in sec t ion [5.31 

ICui et al~l (|2005S) determine a photoabsorption rate 
based on the population of H 2 in the J = and J = 
4 states, the latter of which is optically thin and there- 
fore measures the intensity of radiation outside the cloud. 
Solving for this radiation field they derive J v {\ = lOOOA) 
« 2.1 xlO" 3 J„,q (A = 10001), where J V , Q (A = lOOOA) 
is the UV intensity in the solar neighborhood at 1000 
A and is, according to their paper, equal to J u ,q (A = 
1000A) pa 3.2xl0" 20 ergs cm" 2 s" 1 Hz" 1 sr" 1 . They 
do recognize that the radiation field as determined by 
the H2 J states is ^3 orders of magnitude weaker than 
that determined by the C II* technique and comment 
that this is a discrepancy. More importantly however, 
this radiation field is ~3 orders of magnitude below the 
Haardt-Madau background, which sets a lower limit to 
the radiation field. It is difficult to understand at first 
why the molecular hydrogen excitation is not consistent 
with this minimum radiation field. A solution exists how- 
ever, if the bulk of the H2 gas resides in the narrow C I 
velocity component, i.e. component 2, and hence, has 
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been missed because of the effects of hidden saturation. 
I Car swell et al.l (|2010f ) have completed an analysis of this 
narrow component and include a detailed interpretation 
of the H2 data that is consistent with both the Haardt- 
Madau background and the local stellar held derived by 

cir. 

When we relax the constraint on the radiation held as 
determined by the C II* technique and allow the held to 
vary, as already discussed in section 13.3. 1[ we can con- 
strain the radiation field to j v total & 1 < 1 xlO -19 ergs 
cm" 2 s" 1 Hz" 1 si-" 1 (la). This limit on J v total , c 1 is ac- 
tually quite strong considering that the Haardt-Madau 
background is only ~1 order of magnitude lower than 
this value. Note that the lower limit on the allowed 
j^totai.c i j g fj xec j \yy J v Bkd _ It is apparent from Fig- 
ure [16] that the density is limited to the range, 11 < n(H 
I) < 44 cm" 3 while the temperature is constrained to be 
79 K < T < 794 K, and the pressure 3.5 < log (P/k) 
< 4.04 cm -3 K. We stress that these limits are derived 
independent of the results of the CII* technique, using 
only the CI fine structure data and the assumption of 
ionization equilibrium for a range of possible radiation 
fields. The 2a results provide the less restrictive results 
jjotal < g 6 x iQ-18 ergs cm" 2 s" 1 Hz" 1 sr" 1 , n(H I) 
> 10 cm- 3 , and T < 794 K. 

4.4. J2100-00 

While insufficiently good quality data hinder a full 
analysis of DLA 2100—00, we attempted an analysis of 
component 1, by treating the upper limit on N(CI**) 
as a detection. The measured N(CII*) produces a 
jjocalfiW = 17 x iQ-19 ergs cm -2 s -i Hz -i sr -i ; 

also an upper limit. The relatively large errors on the 
fine structure column densities do not produce very re- 
strictive results. However, the VPFIT derived Doppler 
parameter of this component, b = 0.2 ± 0.3 km s" 1 , 
if real, places a strict upper limit on gas temperature 
of Tthermai ;$ 29 K, which is consistent with the results 
from the CI fine structure analysis. Unfortunately, the 
weakness of the C I line coupled with the low S/N of the 
data, makes a curve of growth analysis inconclusive. For 
the purposes of this paper however, it is the C I column 
densities that are important, and while they typically de- 
pend on the choice of Doppler parameter, we can obtain 
a sort of lower limit to the column density by hxing a rel- 
atively large Doppler parameter, i.e. b = 7 km s" 1 . We 
measure N(CI) = 12.12 cm" 2 , which is well-within the 
errors quoted in Table [2j and see that in this case, the 
C I column densities have little dependence on changes 
in b. 

4.5. Q2231-00 

While there are two C I components contained in the 
DLA at z a i„ s ~ 2.066, here we will focus only on com- 
ponent 2, which possess measurable CI fine structure 
lines. While not formally an upper limit, the detection 
of N(CI**) is relatively weak and has large error. For 
the purpose of this analysis we will treat this value as a 
detection. 

We estimate the jJocai,cn f ronl the measured 
N(C II*) with some conhdence, even though it falls within 
the Lyman-a forest, because the velocity profile of CIT 
closely traces that of other low ions as expected. Using 



this jj°cal,cir ~ 24.7 xlO" 19 ergs cm" 2 s" 1 Hz" 1 sr" 1 , 
we hnd 199 cm" 3 < n < 1513 cm" 3 , and 13 < T < 25 
K, see Table [3 Relaxing the constraint of j v local fiii* 
given by the C II* technique results in a lower limit on 
the density of n(H I) > 3 cm" 3 (2a). 

4.6. J2340-00 

DLA 2340—00 is a relatively complex system requiring 
nine CI velocity components. Given the strong blend- 
ing of transitions and the resultant uncertainty on upper 
limits, we do not analyze those components for which 
N(CI**) is not detected, namely, components 1 and 7. 
In addition, the column densities of component 6 have 
extremely large errors, making analysis pointless, as well 
as the unphysical condition of N(C I*) > N(C I). 

This DLA contains a relatively high total column den- 
sity of molecular hydrogen, logN(H2) = 18.20 cm" 2 , or 
f = 0.014. The neutral hydrogen column density, at log 
N(H I) = 20.35 ± 0.15 cm" 2 , is close to the threshold 
defining a DLA (log N(H I) = 20.30 cm" 2 ), w hich sup- 
ports the findings of iNoterdaeme et all (|2008| ) that the 
probability of detecting H 2 does not strongly depend on 
N(H I). Furthermore, the low- ion velocity prohle is large, 
with Av = 104 km s" 1 , and the cooling rate, as deter- 
mined by CII* over the entire profile, at log £ c — —26.15 
ergs s _1 H _1 , is among the highest of DLAs. Because of 
the complex nature of the CI, the H2, and the low-ion 
prohles, all requiring multiple components, and because 
of heavy blending and saturation in some components of 
the low-ions, we were not able to obtain unique fits to all 
low ion components using VPFIT. Instead, we analyzed 
the CI data using two different model assumptions. 

In the first case, we assume J v total just slightly above 
the background, J u total — 3.9 x 10" 20 ergs cm" 2 s" 1 
Hz" 1 sr" 1 . Because this is a lower limit to the radiation 
field, this analysis provides an upper limit on density (i.e. 
for all else being constant, if we increase the radiation 
field, the density required to collisionally excite the CI 
fine structure levels is decreased). The metallicity, dust- 
to-gas ratio and logoff- are determined by AODM over 
the entire profile with measured results given in Table [H 
The densities and temperatures derived from the C I hne 
structure data are given in Table [SJ 

In an attempt to refine the model, in the second case, 
we used the AODM over three 'super-components' of the 
CII*, S II, Fe II, and other low-ion and resonance line 
transitions. The choice of super-components was moti- 
vated by visual inspection of the spectra that reveal three 
large velocity components separated in velocity space (al- 
beit with each component containing smaller substruc- 
ture). We arbitrarily chose v = km s _1 at z a b s = 
2.054151, and defined the AODM super-components as 
follows: super-component (a) from v = —30 to 15 km s" 1 
, super-component (b) from v = 15 to 70 km s" 1 , and 
super-component (c) from v = 70 to 120 km s" 1 , see Fig- 
ure 1201 Super-components a, b, and c coincide with C I 
components (1, 2), (3, 4, 5, 6, 7) and (8, 9) respectively. 
We then applied the metallicity, dust-to-gas ratio, and 
N(CII*) measurement derived from the AODM 'super- 
component' to each associated CI component. Table [S] 
contains a summary of this analysis. Looking at each 
super-component, we hnd the following: 

•Super — component (a) : We did not perform the 
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C II* analysis on super-component (a) because we mea- 
sure a super-solar [Fe/S] and [Ni/S]. The absence of de- 
pletion detected for Fe and Ni implies that we cannot 
calculate the dust-to-gas ratio. If real, this super-solar 
(or nearly solar) value of Fe II would require a different 
and special star formation history. Given that we cannot 
calculate a dust-to-gas ratio or j v local > cn ; anc j the fact 
that the C I components 1 and 2 associated with super- 
component (a) contain the smallest amount C I, we did 
not perform further analysis. 

• Super — component (b) : Super-component (b), 
which includes CI components 3, 4, 5, 6 and 7, contains 
the bulk of the gas. The C II* technique results in 
j^totaLCir = 5.27xlCT 18 ergs cm" 2 s" 1 Hz" 1 sr" 1 . 
The combination of the large J v total and the low Qjf 
are not compatible with the measurements of C I fine 
structure in components 3, 4, or 5 (6 and 7 are ruled 
out by upper limits). 

•Super — component (c) : The C II* technique applied 
to super-component (c) , which covers C I components 8 
and 9, results in a J v total = 1.1 xl0" 18 ergs cm" 2 s" 1 
Hz" 1 sr -1 . A reasonable range of results for components 
8 and 9 are summarized in Table [5j 

If we relax the constraint of J v t ° tal > c 11 as measured 
by the C IT technique, we can determine limits for each 
component for a range of possible radiation fields. Re- 
sults are given in Table [6l In all cases we can put a 
lower limit on the density of n > 7 cm -3 . In all com- 
ponents except for component 2, the temperature can 
be constrained to be T < 800 K and the J v total ^i < 
27 xl0 _19 ergs cm" 2 s" 1 Hz -1 sr" 1 , is 2 orders of mag- 
nitude larger than the Haardt-Madau background and 
consistent with the J v total derived by the C II* technique. 

4.6.1. On the possible ionization of DLA 2340—00 

Given the potential presence of Fe III A1122 and the 
fact that the column density of DLA 2340-00, at log N(H 
I) = 20.35 cm -2 , is just above the DLA threshold of 2 x 
10 20 cm -2 , we considered the possibility of partial ion- 
ization of the gas. We used the AODM technique to mea- 
sure Fe III/Fe II in the same three 'super-components' 
discussed above, covering CI components 1 and 2, com- 
ponents 3, 4, 5, 6 and 7 and components 8 and 9 respec- 
tively, see Figure Because the Fe III profile does not 
trace that of the low ions (see Figure l2Tfl) . we cannot rule 
out the possibility of blending with the Lyman a forest 
and therefore, we treat our Fe III measurements as upper 
limits that ultimately give no information about possible 
ionization. 

Instead, we can rule out a high ionization fac- 
tor based on the [A r/Sl measurement. Specifically, 
iProchaska et al.l (|2002f ) state that photoionzation models 
in which [Ar I / S II] > -0.2 dex require that x < 0.1, 
in other words, require gas that is > 90% neutral. Using 
the Ar I AA 1048, 1066 transitions, excluding blending 
in the super-component 3 of the A1066 transition (see 
Figure |2"0")) . we find [Ar/S] > —0.2 in all components (see 
Table O, indicating that x < 0.1 and that the gas is 
>90% neutral. 

4.6.2. Molecular Hydrogen tn DLA 2340- 00 

The total column density of molecular hydrogen, 
logN(H 2 ) = 18.20 cm" 2 , where f = 0.014, is large rel- 



ative to most H2 -bearing DLAs, where f is typically f ~ 
10~ 5 . We have analyzed the H 2 using VPFIT. To allow 
for the best fit we have let the H2 component redshifts 
and b values vary independently of the C I and low-ion 
components. Wc find that wc require 6 H2 components to 
achieve the best fit. In redshift space, these components 
lie remarkably close to the CI components 1, 2, 4, 6, 8 
and 9, (Aw = +1.4, +0.6, -2.8, +0.3, -0.4, +0.4 km s" 1 
, respectively) and therefore we use this notation to refer 
to the H2 components. Details of the H2 measurements 
are given in Table [TU] and example spectra in Figure 1211 
In Figure [H] we plot the excitation diagrams for each H2 
component and list the excitation temperature as deter- 
mined by the J=0 and J=l states. Additionally, we use 
the population of the J=4 state to determine the radia- 
tion field as described in Appendix 3. Details are given 
in Table QT] 

While the radiation fields derived from the H2 data are 
in general consistent with the C I constraints, it is inter- 
esting to note that the densities derived from the H2 data 
alone, see Appendix 1 for details, while consistent with 
the C I limits, tend to be significantly higher than that 
required for C I. We summarize results for each compo- 
nent for which we could make comparisons between the 
different techniques: 

• Components 1& 2 : Excitation temperatures are in 
good agreement with constraints placed by the C I data, 
however the H2 analysis was not completed because of 
the super-solar [Fe/S] and [Ni/S]. 

• Component 4 : The H2 derived T = 276 K is consis- 
tent with the 2a C I range of T = 40 - 794 K. The density 
derived from H2, n(H I) ^1600 cm" 3 is not compatible 
with the C I limits. Finally, the H 2 derived JJ"*al,H 3 = 
2.25 x 10" 19 ergs cm" 2 s" 1 Hz" 1 sr" 1 is within range 
allowed by C I and much lower than that predicted by 
GIF ( jJotai,cir = 52 7 x 10 -i9 ergs cm -2 g -i Hz -i 

sr- 1 ). 

• Component 6 : Because we only obtained upper lim- 
its on the C I** state of C I component 6 we can not make 
direct comparisons in this case. We can however, com- 
pare with CI component 5. In this case, the H2 derived 
T = 587 K is consistent with the 2a C I limit T < 1259 K 
(or T < 794 depending on J„ ), and the H2 derived n(H 
I) = 10, 509 cm -3 is also consistent with the lower limits 
placed by C I data. The radiation field derived from the 
H 2 , jJ° tal ' H i = 18.8 x 10" 19 ergs cm" 2 s" 1 Hz" 1 sr" 1 , 
is consistent with the limits placed by C I, however these 
are both much lower than the j^ total < CH predicted by 
the C II* technique. 

• Component 8 : The H2 derived T = 475 K is incon- 
sistent with the 2a CI limit, T < 158 K. The H 2 derived 
density, n(H I) = 3595 cm" 3 , is also incompatible with 
the range allowed by C I (n(H I) = 12 - 209 cm" 3 ) . How- 
ever, the H2 derived radiation field, j v total - R i < 5 13 x 
10~ 19 ergs cm" 2 s" 1 Hz" 1 sr" 1 is compatible with the 
limits set by C I, and similar to that predicted by the 
CII* technique, J v toi < c n * = 11 x 10" 19 ergs cm" 2 
s _1 Hz" 1 sr" 1 . 

• Component 9 : The H2 derived T = 151 K is consis- 
tent with the 2a CI limit, T < 398 K. The H 2 derived 
density, n(H I) = 377 cm -3 , is also nearly compatible 
with the range allowed by CI (n(H I) = 28 - 363 cm" 3 ). 
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However, the H 2 derived radiation field, JJ° tal > H 2 < 0.48 
x 10~ 19 ergs cm~ 2 s _1 Hz -1 sr -1 , while compatible with 
the limits set by C I, is not compatible with that pre- 
dicted by the CIF technique, j u total ^ir = n x 1Q -i9 
ergs cm -2 s _1 Hz -1 sr _1 . 

A striking difference is seen in the densities de- 
rived from the CI and H2 data. A similar differ- 
ence has been reporte d in sightlines toward s Q001 3— 004 
and Q1232+ 082 (IHirashita fe Ferrara I 120051) an d 
HE0515-00 (jReimers et al.l l2003t iQuast et al. I [2001 . 
They argue that if the density is this high, i.e. equal 
to the critical density, then the states should be ther- 
malized and there would be no difference in the excita- 
tion temperature as determined by the J = 1 and J = 
2 states. Given t hat the observed T02 is alw ays much 
greater than T01, IHirashita fc Ferrara I (120051) p ropose, 
following the suggestion of iReimers et al.l (|2003T ) , that a 
potential explanation may be that the H2 formation rate 
O^dust) may be larger than estimated, which would result 
from, for example, a smaller than estimated grain size. 
Therefore, if the formation rate is higher, you require a 
smaller density than that derived. On the other hand, 
for each component of DLA 2340—00 studied here, the 
T02 is either less than or approximately the same as T01 , 
implying that the critical density cannot be ruled out. 

5. DISCUSSION 

We have used C I fine structure absorption in high res- 
olution, high signal-to-noise data to study the physical 
conditions in DLAs at high redshift. Our work differs 
from previous studies of C I fine structure absorption be- 
cause we did not assume a gas temperature in order to 
derive the density. Rather, we assume ionization equi- 
librium, which in conjunction with the C I fine structure 
data, allows us to constrain both the temperature and the 
density of the cloud. In addition, we use the C II* tech- 
nique to infer the local radiation field due to stars and 
include its contribution to the CI fine structure excita- 
tion, thus providing a complete and fully self-consistent 
model of the gas. In most cases, the CI fine structure 
excitation is consistent with the j v local > c 11 ; derived in- 
dependently by the C II* technique. 

5.1. Assessment of Systematic Errors 

To draw meaningful comparisons between the CI re- 
sults and those of the C II* technique we must first ana- 
lyze our possible systematic errors. The primary source 
of error in the CI analysis stems from two sources - 1) 
measurement error of the fine structure column densi- 
ties, including possible errors in the oscillator strengths 
(we report results for 2a errors, which should encompass 
these errors 12 ), and 2) the assumption that in solar units 
the carbon abundance is equal to that of Si II (or S II) 
—0.2 dex in case of minimal depletion. The source of 
systematic errors in the CII* technique are more diffi- 
cult to assess (see WGP03 for a complete discussion). 
We explore the possible effects of these errors in the fol- 
lowing section and demonstrate that, even considering 

12 For example, tests run using Jenkins 2006 f- values for C I 
produced a Doppler parameter ~ 15% larger and column density 
~0.3 dex smaller than the Morton 2003 values in the case of the 
curve of growth analysis of the narrow, cold component in DLA 
0812+32 [IJorgenson etaIll200SD . 



these errors, the CI data appear to be probing gas of 
higher densities and pressures - likely small knots of gas 
within the larger DLA galaxies - than that probed by 
the 'global' C II* technique. 

5.1.1. Comparison with CIF technique model 

While it is difficult to assess the systematic errors in- 
volved in the C II* technique, some of the potential un- 
certainties have been removed since the work of WGP03. 
For example, the SMC dust model is now assumed to be 
correct because of the non-detection of the 2 175 A dust 
feature that would have indicated Galactic dust, while 
the redde ning curve resemble s SMC rather than Galactic 
orLMC l|Vladilo et all2 008). Therefore, we focus on the 
two largest potential uncertainties: 1) the assumption of 
the equilibrium pressure, P eg , and 2) the minimal ver- 
sus maximal, depletion model (see discussion in section 
2.3 of WGP03). The standard CII* technique involves 
solving the equations of thermal a nd ionization e quilib - 
rium as described in WPG03 and iWolfire et ail (fl995l) . 
A unique solution is determined by assuming that the 
equilibrium pressure is equal to the geometric mean be- 
tween V max and P min P P — CP P W 2 

lwLLn r aliu r , r e q — r geo — \± m%n r max ) 5 

where P m in and P mM are the minimum and maximum 
pressures of the function P(n) where n is gas density. 
This results in two stable solutions for a given star for- 
mation rate, one WNM and one CNM, and is the basis of 
the two-phase model. However, a two-phase medium can 
achieve equilibrium with a pressure ranging from P m in to 
Pmax, and therefore, the assumption of P eq equal to the 
geometric mean, while reasonable (see Wolfire et al. 2003 
and WPG03 discussion), is still an unproven assumption. 
Following WGP03 we attempt to gain a sense of the pos- 
sible systematic errors by allowing P e? to vary between 
Pmin and P m ax- In Figure [23"1 we show the standard two- 
phase diagram, plotting in (a), log(P/k) versus density 
for various star formation rates per unit area (which is 
proportional to J„) which are constant along each P(n) 
curve, and which increase from bottom to top, and in 
(b) log(^ c ) versus density for £ c equilibrium solutions for 
those same star formation rates. The green dashed line 
in each plot indicates the P(n) solution for heating by 
background radiation alone (i.e. logSsFR — — oo). The 
black horizontal line in (b) denotes the observed cooling 
rate of DLA 1331+17, log(4) < -27.14 ergs s" 1 H _1 . 
Three vertical, black, dotted lines illustrate the location 
of the CNM stable points associated with, from left to 

right, Pmin: P geo (P minP max ) ^ 3 and P max for the 

case of background radiation alone. The stable (£ c , n) 
pairs for the grid of star formation rates are denoted by 
the three red lines in (b), where the three different pres- 
sure assumptions have been made - from left to right 
they are: dashed = P m m, dot-dashed = P geo and dotted 
= Pmax- Although not relevant for our current discus- 
sion, it is seen that P m in requires a higher star formation 
rate and lower density to achieve an equilibrium solution 
with a cooling rate equal to that observed. On the oppo- 
site extreme, P m ax requires a higher density and lower 
star formation rate. 

For the purposes of this paper, we are interested in 
the range of densities and temperatures that result from 
these different model assumption inputs to the C II* tech- 
nique. In Figure [24] we plot the resultant log(P/k) ver- 
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sus density for our example case of DLA 1331+17. For 
the minimal depletion model, the case of P eq — P geo = 
(PminPmax) 1 is denoted by the asterisk. Lines con- 
nect to the P m in and P ma i solutions, providing a sense 
for the potential systematic "error" involved in the as- 
sumption of where the equilibrium pressure resides. We 
have also plotted the results of the maximal depletion 
model, denoted by the diamond, with dashed lines con- 
necting to the associated P m in and Pmax solutions. To 
correctly compare this range of solutions to that of the 
C I data we must re- model the C I theoretical curves in 
each case because the change in the assumed P eq results 
in a change of star formation rate, or J v local , which is 
an input for the CI theoretical curves. We plot the 2a 
results of the CI analysis for J u local spanning that deter- 
mined by P m in and P m ax- We summarize the results of 
these two techniques for DLA 1331+17 in Table [T2l 

In Figure[25]we summarize the results of our C I sample 
and compare them to the C IP-derived and H2-derived 
models of the same systems. We plot the best- fit C I solu- 
tions, as determined by x 2 minimization, as circles, while 
the shaded regions denote the 2tr C I solutions. The C II* 
solutions are represented as diamonds. It is seen that the 
densities determined by the C II* technique are systemat- 
ically lower than those determined by the C I data. As a 
result, the overall pressures are lower. The temperatures 
vary depending on model assumptions, but are in general 
agreement or higher than the CI results (see Table [T2")) . 
We conclude that the C I gas is tracing a denser region of 
the DLA than that traced by the global CII* technique: 
essentially, the C I resides primarily in small, overdense 
knots. However, one problem with this picture is under- 
standing how these two phases, the low pressure ambient 
medium and the higher pressure smaller 'clumps' of C I 
-bearing gas, remain in pressure equilibrium. We return 
to this issue in a following paper. 

These comparisons also reveal a discrepancy between 
the ^fj- derived by the C II* technique analysis and the 

observed ^jj-- In general, the observed ^tj- ratio is ap- 
proximately an order of magnitude smaller than that pre- 
dicted by the CII* technique model (see Table [T2")) . Wc 
avoid a detailed discussion here and refer the interested 
reader to WGP03, section 5.1, for a detailed discussion 
of the model inputs that affect ^fj-- We only briefly 
mention that we have tested the results with cosmic rays 
turned on and off, and there is not a large effect on the 
results for the star formation rates we are considering. 
However, this is not the case for large star formation 
rates, J v total ^10 -18 ergs cm~ 2 s _1 Hz -1 sr _1 , for which 
the effect of cosmic rays becomes more important. As- 
suming that the C IT technique model is correct, we can 
understand this difference in terms of the new model in 
which the bulk of the measured C I is localized in small, 
dense clumps, « 1 — 10 pc, relative to the larger DLA 
(for estimates of the cloud size, see last column in Ta- 
ble [5]). In this case, while the cloud is still optically 
thin and feels the same radiation field as the surround- 
ing medium, the increased density, and hence, increased 
electron density, n e , of the CI cloud work to lower 

because is inversely proportional to density for a 
fixed radiation field. In a sense, this conflict with the 
predictions of the 'global' CII* technique model, is an 



expected result of the overdense-C I -region model. 

5.1.2. Direct Estimate of Error on jJ ocal 

While there are many assumptions and uncertainties in 
the C II* technique, we can make a direct estimate of the 
possible error on J u local derived from the C II* technique 
by usi ng the z abs =1.9 DLA 2206-19. In IWolfe et all 
(2004), the FUV radiation field infer red from an image o f 
a galaxy associated with this DLA (jWarren et al.ll200l[) 
is compared with the radiation field inferred from the 
C II* technique and they are found to agree to within ~ 
50%. Therefore, we will assume that the error on J v local 
is - 50% . 

5.1.3. Distribution of C II 

Throughout this work we assume that the measured 
N(CII) can be directly applied to each CI cloud. How- 
ever, assuming that there is no hidden saturation - which 
would actually increase the N(C II) - this is actually an 
upper limit to the amount of N(C II) associated with each 
C I cloud. The distribution of C II along the line of sight 
could be clumpy such that only a fraction of the mea- 
sured N(C II) is associated with a given C I cloud. With 
respect to the analysis done in this paper, the effect of 
decreasing the amount of CII associated with the CI 
cloud causes an increase in the resultant volume density 
and a decrease in the temperature such that the pressure 
remains approximately constant. The increase in volume 
density with decreasing fraction of CII results in a de- 
crease of the derived cloud size as shown in Figure 1261 

5.2. Relation to 'high-cool' DLA population 

With the exception of DLA 1331+17, all of the CI - 
bearing objects not only contain C II* (as compared with 
^50% of the general DLA population) but also have cool- 
ing rates, £ c , th at place them firm ly in the 'high-cool' 
range defined by IWolfe et al.l (|2008l ) (median 'high-cool' 
log£ c = —26.6 ergs s _1 H^ 1 ). A literature search re- 
veals that all previously published DLAs with positive 
detections of CI fall into the 'high-cool' population as 
wel l. While it is not cl ear how the bimodality discovered 
by IWolfe et al.l (|2008f ) is related to C I, this correlation 
- that almost all CI -bearing DLAs are also high-cool 
DLAs - could simply be a result of the higher metallici- 
ties and dust-to-gas ratios that are required to form and 
sustain, through dust shielding of UV radiation, mea- 
surable amounts of H2 and C I. This trend, of higher 
metallicity DLAs being more likely to c ontain measur- 
able H 2 , has been observed previously by pctitica n et al.1 
(120061)' Additiona lly, the 'high-cool' DLAs, shown by 
IWolfe et all (|2008f ) to consist of primarily CNM, simply 
might be an environment more conducive to the presence 
of H 2 and C I. 

5.3. Comparison with the local ISM 

In this section we compare our C I-bearing clouds to 
the ISM of the local Universe, namely, the Milky Way 
and the Small and Large Magellanic Clouds (the SMC 
and LMC respectively). To draw these comparisons, we 
first determine the median n(H I), T, and P of our CI 
sample. Of course, this median is dependent upon which 
models we choose to include. Here, we include the models 



Using C I fine structure to probe physical conditions in DLAs 



15 



for each DLA that are likely to be the most physically re- 
alistic, i.e. case (a), the \ow-£ c , minimal depletion model 
for DLA 0812+32, z abs =2M, and the 'global' models 
for both DLA 0812+32, z abs = 2.62 and DLA 2340-00. 
We chose the 'global' models as the most likely to be 
correct because of the inherent uncertainty in determin- 
ing the distribution of both the N(H I) and the met- 
als amongst the velocity components necessary for the 
component by component analysis. A summary of these 
models is shown in Figure 1251 The resulting median val- 
ues for this sample are: <n(H I)> = 69 cm" 3 , <T> = 
50 K, and <log(P/k)> = 3.86 cm" 3 K , with standard 
deviations, o n {H i) — 134 cm" 3 , <jt — 52 K, and <Ji og (p/k) 
= 3.68 cm- 3 K . 

In the Milky Way, Uenkins fc Tripp I (|2001l ) used very 
high resolution (R = 200,000) STIS data to analyze C I 
fine structure populations and find a median log(P/k) = 
3.35 cm" 3 K . This median pressure is similar to the pres- 
sures we derive in high redshift C I -bearing DLAs, where 
the median pressures derived are typically log(P/k) = 3 
— 4 cm -3 K . While this would seem to indicate that 
physical conditions similar to the Milky Way exist in 
high-z DLAs, we point out that in the case of the Milky 
Way, the high pressure is driven by the much higher heat- 
ing rate (log€ c ~— 25 ergs s" 1 H" 1 ), a result of the higher 
dust-to-gas ratio (typically 30 times that of DLAs). In 
contrast, the heating rates in high redshift DLAs are gen- 
erally 1—2 orders of magnitude smaller. Hence, the pres- 
sures derived in the C I -bearing clouds, while similar to 
those observed in the Milky Way, are not the result of the 
same physical condition s as th ose observed in the Milky 
Way. Uenkins fc Tripp 1 ()200lD also observe a small pro- 
portion of the gas in many sightlines to be at very high 
pressures, P/k > 10 5 cm -3 K , which they speculate are 
caused by converging flows in a turbulent medium or in 
turbulent boundary layers: such pressures have not been 
detected in DLAs. 

It is perhaps more meaningful to compare DLAs with 
the Large and Small Magellanic Clouds that, like DLAs, 
are known to have sub-Milky Way du st-to-gas ratios and 
metallicities. iTumlinson et al.l (|2002l ) performed a FUSE 
survey of H2 along 70 sightlines to the Small and Large 
Magellanic Clouds. For all sightlines with logN(H2) > 
16.5 cm" 2 they find < T 01 > = 82 ± 21 K, whereas for all 
sightlines (including those at lower columns), they find 
< Tqi > = 115 K. This can be compared with the Galac- 
tic average, T = 77 ± 17 K (|Savage et alj|1977l) . The 
temperatures found by the CI and H 2 data presented 
here are in broad agreement with these values. A notice- 
able exception to this agreement are the temperatures 
derived from the H2 rotational states in DLA 2340—00. 
They are generally higher (T ex 01 w 150 - 600 K) than 
those found in the SMC/LMC and in better agreement 
with the mean ki netic temperature of th e gas, T = 153 
± 78 K, found by iSrianand et al. I ([2005) in a sample of 
H2 -bearing DLAs at high-z. 

6. CONCLUSIONS 

The goal of this paper is to present new detections of 
neutral carbon in high redshift DLAs and to present a 
new method for analyzing the C I fine structure lines. As 
done by several other authors, we utilize C I fine struc- 
ture lines to determine densities, however, instead of as- 



suming a temperature, our work constrained the allowed 
density and temperature combinations using only the col- 
umn density of C I in the fine structure states and the 
assumption of ionization equilibrium. In a second paper, 
we will incorporate these physical conditions into a gen- 
eral model for C I-bearing DLAs. Our major conclusions 
are as follows: 

1. The steady state analysis of CI fine structure pop- 
ulations along with the assumption of ionization 
equilibrium provides realistic constraints on both 
the volume density and temperature of high red- 
shift DLAs. The C I data are in general consistent 
with the radiation fields, jJ ocal - c 11 ; derived from 
the C II* technique and provide further evidence of 
the presence of CNM in high redshift DLAs. 

2. The densities and pressures of the CI- bearing gas 
are systematically higher than those of the 'global' 
DLA predicted by the C II* technique. We propose 
two physical scenarios that could be consistent with 
the data presented here. First, the CI could be 
tracing overdensities that are created by shocks, 
hence the CI exists in the post-shock cool gas. 
However, it seems likely that the post-shock gas 
would have a systemic velocity offset from the pre- 
shock gas. This is not observed. In fact, the general 
good agreement between the velocity centroids of 
the C I and other resonance lines argues against the 
shock idea. A second scenario is that the C I exists 
in higher density, higher-pressure edge of a pho- 
todissociation region, i.e. the edge of a molecular 
cloud. While the high optical depth through a clas- 
sical molecular cloud would obscure a background 
quasar, the photodissociation region, or edge of the 
molecular, cloud could be optically thin enough to 
allow transmission of the background quasar light 
and would be consistent with the gas physics de- 
termined by the C I fine structure lines. 

3. As noted by ISrianand et al. I (|2005l ). all CI - 
bearing DLAs also contain CII* absorption. We 
find that, with only one exception, all CI objects 
for which C II* coverage is available (5 presented in 
this paper, 1 yet to be published, and 7 from the 
literature) contain strong C II* absorption, placing 
them in the category of 'high-cool' DLAs. This 
could be simply a consequence of the fact that C I 
-bearing DLAs generally host larger fractions of H2 
whose formation is encouraged by the higher than 
average metallicities and dust-to-gas ratios, consis- 
tent with the 'high-cool' population of DLAs. 

4. High resolution studies of neutral carbon lines re- 
veal narrow, sub-lkm s _1 , cold and unresolved 
components. These components likely contain rel- 
atively large amounts of gas and are most likely 
cold, dense knots, perhaps photodissociation re- 
gions on the edges of star forming regions. This 
would explain the presence of CI, H2, and larger 
than average dust-to-gas ratios. To date, two such 
components have been published; in DLA 0812+32 
with a temperature of T < 78 K (|Jorgenson et al.l 
[20091 ) and in DLA 1331+17 with T < 218 K 
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([Carswell et al.l I2010D and we presented an addi- 
tional candidate in this paper. Such clouds may 
exist in all C I systems. Their non-detection in 
other C I systems does not rule out their existence 
due to the difficulty in detecting such small equiv- 
alent widths that are likely blended with other ve- 
locity components. It is possible that these narrow 
components are ubiquitous and contain significant 
amounts of gas that has been previously missed in 
lower resolution studies. 
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7. APPENDIX 1: MOLECULAR HYDROGEN 

The purpose of this Appendix is to present an outline of 
the molecular hydrogen analysis performed in this work. 
While the focus of this work was the analysis of neutral 
carbon, DLAs that contain neutral carbon are likely to 
also contain detectable H2. This is because neutral car- 
bon and molecular hydrogen are photoionized and pho- 
todissociated respectively, by photons of the same ener- 
gies and therefore, they are usually found together. This 
is the case in several of the DLAs presented in this work, 
i.e. DLA 1331+17, DLA 0812+32, and DLA 2340-00. 
For DLA 2231-00 and the low-z DLA 0812+32, we did 
not have coverage of the H2 region and therefore cannot 
determine anything about the presence of H2. Here we 
present an outline of the DLA H2 analysis frequently 
perform ed in the literature, m os t recently by works 
such a s ILevshakoy et al. 1 (I2002D. iHirashita fe Ferrara I 



(l2005lh iCui et al. I (l2005h. iNoterdaeme et al. I (|2007aD. 
and Noterda eme et al. Trt2qo7bn and based upon work 
bv ISpitzer fe Zweibell (|1974D . Uural (j!975aft and I Jura I 



(|1975bD . Essentially, the measurement of H2 in the differ 
ent rotational J states allows for an independent estima- 
tion of the physical properties of the cloud such as den- 
sity, temperature and the incident radiation field. While 
the H2 analysis arguably contains several uncertainties 
and assumptions, it is nonetheless interesting to compare 
the results from these two independent methods, the H2 
analysis we present below and the C I fine structure anal- 
ysis presented in this work. 

First, we can estimate the kinetic temperature of the 
cloud using the column densities of H2 in the J=0 and 
higher J rotational states. Ac cording to the Bo l tzman n 
distribution (see equation 8 in ILevshakoy et al.l ((2002)): 



NjJ) _ g{J)_ _ B V J(J + 1) 
N(0) ~ 5 (0) 6 T ex 



(8) 



where = 85.36 K for the vibrational ground state and 
g(J) is the degeneracy of level J, given by, for level J = x, 
g x = (23 x + 1) (21^ + 1) where I = for even x and I 



1 for odd x (For states J = 0-5, g = [1, 9, 5, 21, 9, 33]). 
The excitation diagram is typically plotted as log (Nj / 
gj) versus the relative energy between the level J and J 
= 0. The excitation temperature, defined in equation 
is inversely proportional to the negative slope of the line 
connecting the excitation diagram points, i.e. 
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The typical assumption is that the kinetic temperature 
of the cloud can be estimated by the excitation temper- 
ature derived from the population of H2 in the states J 
= and J = 1, assuming that the J = 1 level is thermal- 
ized. Because the critical density of the low J states is 
relatively small, the population of the low J levels is gen- 
erally dominated by collisional excitation and therefore 
it reflects the kinetic temperature of the gas when in lo- 
cal thermodynamical equilibrium. The higher J states (J 
> 2) are typically characterized by a higher T ex , or flat- 
ter slope, which is explained by population mechanisms 
other than collisions. 

The higher T ex derived from the population of H2 
in the higher rotational J state s was originally unex- 
pected (jSpitzer fe Zweibelll974f ). After observations by 
Spitzer and Chochran in the 1970s in which they ob- 
served the hig h excitation temperatures derived from the 
high J states, ISpitzer fe Zweibell (|1974fl proposed meth- 
ods other than collisions that could populate the higher 
J states. They proposed two methods other than colli- 
sions: 1) the cascade down from upper vibrational levels 
following the absorption and reemission of Lyman and 
Werner band photons -i.e. the J = molecule is ex- 
cited to a higher vibrational state and then de-excites to 
a higher J state in the ground vibrational state (rather 
than staying in the original J=0 or J=l state), or 2) the 
direct formation of H2 in a higher state (i.e. the H2 pops 
off the dust particle in an excited state and cascades to 
say J = 4). Since these processes are believed to domi- 
nate collisions, the typical practice is to neglect collisions 
and to consider only these two processes, as we show in 
the following analysis. 

Assuming steady state and neglecting collisions, we use 
the two afore mentioned populating mechanisms and as- 
sume depopulation by spontaneous emission (true until 
the density is above 10 4 ), to write the following steady 
state equations for state J = 4 and J = 5 respec tively 
(these are equations 2a and 2b from I Jura I (|1975b| )). 

p 4 oPon{H 2 , J = 0)+0.19Rn(HI)n(H) = A 4 2 n(H 2 , J = 4) 

(10) 

and 



p 5 ,oPin(H 2 , J = l)+0.URn(HI)n(H) = A 5 , 3 n(H 2 , J = 5) 

where p is the pumping coefficient or pumping efficiency 
(or by Jura, the redistribution probability) into the J=4 
and J=5 levels from the J=0 and J=l levels respectively, 
/3 is the photoabsorption rate, R is the H 2 molecule for- 
mation rate (on dust grains, also known as Rdust), n(H 
I) is the neutral hydrogen number density, n(H) wn(H I) 
+ 2n(H2), and A is the spontaneous transition probabil- 
ities. In other words, the first term is the UV excitation 
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and decay to higher J state, while term 2 represents the 
direct formation in the higher J state. The values of the 
constants are as follows: A4.2 = 2.75 x 10 _9 s _1 , A 5i3 = 
9.9 x lO-V - 1 (Spitzer 1978), p 4 , = 0.26, p 5 .i = 0.12 
(jJura 111975b! ) . Solving for /3 will allow us to estimate the 
UV field incident on the cloud, while /3 together with R 
will allow us to estimate the neutral hydrogen density. 

In order to solve for j3 and to get rid of the N(H I) 
dependence in equation 1101 we take advantage of the 
assumption of equilibrium between the formation and the 
destruction of H2 (which is reasonable because the time 
scales of H2 formation and destruction are well below a 
dynamical time), as follows, 



Rn(HI)n(H) = R dlss n{H 2 ) 



(12) 



If we substitute equation [T^Jinto equation 1101 and make 
the common assumption t hat 11% of photoabsorption 
leads to photodissociation ()Juralll974bT) , 



R dlss = 0.11/3 



(13) 



we obtain equilibrium equations that are independent of 
the neutral hydrogen column density Ajji, 



correction for the helium content, and id is the cross- 
section of a grain, a d = 7r a 2 . ad = tuv is the optical 
depth in dust and is expressed by iHirashita fc Ferrara I 
(pfiflh as, 



_ 4.2N H m H D 

T uv = jts 



0.879( 



0.1/im 



2g cm 



10 21 cm- 



(18) 



where a is the radius of a grain, S is the grain ma- 
terial density, and D is the dust-to-gas mass ratio. 
IHirashita fc Ferrara I (p005l ) assume the Galactic (Milky 
Way) dust-to-gas mass ratio to be D© = 0.01. They 
define the normalized dust-to-gas ratio K = D/Dq. As- 
suming a=0.1 and <5=2, equation 1181 can be written, 



tuv = 0.879k( 



H 



10 21 cm 2 ' 



(19) 



see iCui et af~l (|2005[ ) equation 7. 

The first part of equation [TTl the self-sheidling 
of H2, is taken from an analytic approximation from 
IDraine fc Bertoldi I (fl996l) and is valid for N(H 2 ) > 10 14 
cm -2 . We can therefore rewrite equation [17] as, 



P4,oPo 



N(H 2 ,J = 0) 
N(H 2 ) 



0.021/3 - A A2 



N(H 2 ,J = A) 



N(H 2 



and 



(14) 



N(H 2 ,J = 1) 
N(H 2 ) 



0.049^i = A 53 



N(H 2 ,J = 5) 
N(H 2 ) 



(15) 

Using the measured H2 column densities we can then 
solve for /3, the photoabsorption rate of H2 in each com- 
ponent. Note, that /3q should equal /3%. 

Once we have solved for /3 we can determine the in- 
cident radiation field by using the relation between the 
radiation field and the photodissociation ra te that it in- 
duces on th e molecular hydrogen. Follow ing lAbel et al. I 
(|1997t) and IHirashita fc Ferrara I (|2005l ). 



Rdiss = (4vr)l.l x 10 8 J^ w S shleld s 



(16) 



where Rdiss is the photodissociation rate (= 0.11/3 as 
above), J V LW (LW stands for Lyman- Werner) is the UV 
intensity at \\v = 12.87 eV averaged over the solid angle 
(12.87 eV is the dominant energy at which the photodis- 
sociation happens). This can be compared with the J v 
that is typically calculated in the C II* technique at A 
= 1500A, or 8.27 eV). S s hield accounts for shielding due 
to two effects: 1) dust shielding and 2) self-sheilding. In 
order to solve for J V LW we must determine the effects 
of shieldi ng in equation IT^l W e estim ate the shielding, 
following Hirashita fc Ferrara I (|2005l ). as 



S 



shield 



(• 



N{H 2 ) 



-0.75 e -a d N d 



(17) 



where the first term expresses the self-shielding and the 
exponential term is the shielding due to dust. Nj, the 
column density of dust, is related to the HI column den- 
sity by: (4/3)7r a 3 £N d = 1.4 mjjN/jD, where 1.4 is the 



S 



shield 



N{H 2 ) , 
10 14 cm- 2 ' 



-0.75 



A, 



esp[-0.879*( " _ 2 )] 

(20) 

Therefore, given /3 and S s hieid we can use equation [16] to 
solve for the ambient radiation field, J U L . Note that 
this is the total radiation field, or 3 V LW = J u total . 

We can also use the measurements of N(H2) to esti- 
mate the volume density of H I. We define the molecular 
fraction, fu 2 as follows, 



fn-2 



2n(H 2 ) 



2N(H 2 ) 



n(HI) + 2n(H 2 ) N(HI) + 2N(H 2 ) 



(21) 



where we measure the N(H I) and the N(H2) directly. 
If we substitute [21] into [12] and remember that n(H) 
wn(H I) + 2n(H2), we can solve for the number density 
of hydrogen, 



Rn(HI) = R d , 



n(H 2 ) 



or 



n(HI) + 2n(H 2 ) 

Rdiss Jh 2 



Rdi. 



/h 2 



n{HI) 



R 



(22) 



(23) 



This is temperature dependent however, and we 
will use the detai l ed ex pression for R given by 
IHirashita fc Ferrara I J2005) (note, they call it R^ S f), 



R = 4Txl0- 17 5 d( r)(-^-)- 1 ( -^)(-^-)V 2 ( 



'O.ljum' v 10- 



100A y 



2gcm 3 ' 
(24) 

where S<j(T) is the sticking coefficient of hydro- 
gen atoms onto dust and everything else was de- 
fi ned previously. The stickin g coefficient is g iven by 
C IHollenbach fc McKee I <fl979h : IQmukai I (12000) 
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S d (T) = [1 + 0.04(T + T d ) - 5 + 2 x 10~ 3 T + 8 x 10- 6 T 2 ]- X x 
(1 + exp[7.5 x 10 2 (l/75 - 1/Td)])- 1 

(25) 

where is the dust temperature and is given by 



(26) 

where A is a constant that depends on the optical prop- 
erties of the dust grains. For silicate grains, A = 1.34 x 
10~ 3 cm and fo r carbonaceous grains , A = 3.20 x 10~ 3 
cm. Following Hiras hlFa fc Ferrara I (|2005f ) we assume 
Quv = 1 (Quv is the dimensionless absorption cross- 
section normalized by the geometrical cross-section) , A 
= 3.20 x 10~ 3 cm and a = 0.1 /im, while x> the nor- 
malized radiation field, was calculated previously from 
the H 2 levels ( X = J„ LW / J u LW e where J„ LW Q = 3.2 
x 10 -20 ergs cm' 2 s" 1 Hz' 1 sr" 1 . Therefore, we can de- 
termine R as a function of T. In the present work, we 
assume that T is equal to the excitation temperature as 
derived from the J=0 and J=l H2 states. Finally, we use 
equation [531 to estimate the neutral hydrogen density. 
Note, this method of determining temperature, density 
and radiation field, is independent of the C I fine struc- 
ture data. 
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Fig. 1. — DLA 0812+32 CI velocity structure. Spectral regions covering the five CI multiplets used in the analysis of DLA 0812+32: 
The multiplets are a) 1656A, b) 1560A, c) 1328A, d) 1280A, and e) 1277A. Black is the data, red is our fit. The different fine structure 
transitions are color coded: CI = red, CI* = green, CI** = blue, while the three velocity components are denoted by different linestyles: 
component 1 at v ~ —43km s~ 1 , or z a ^ s = 2.625808 is dashed, component 2 at v = 0km s — 1 , or z a ^ s = 2.6263247 is solid, while component 
3 at v ~ +14km s — 1 , or z a (, s = 2.626491 is dotted. Interloper lines are denoted by dotted black lines. 
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Fig. 3.— DLA 0812+32 

z a6s — 2.06 CI velocity structure. Notation is as in Figure 1. Only the 1656A and 1560A multiplets were used 
in the fit because of the serious blending with the Lyman a forest in lower wavelength multiplets. The single CI component is located at 
v = km s" 1 , with z abs = 2.066780. 
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Fig. 5. — Spectral regions covering the five CI multiplets used in the analysis of DLA 1331+17. Notation is as in Figure 1. This system 
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Fig. 7.— Spectra of DLA J2100-06 covering the 1656A, 1560A, 1328A, 1280A, and 1277A multiplets. Notation as in Figure 1. The best 
fit requires three CI velocity components. The lower S/N in 1656A multiplet (panel a) is due to its proximity to the order gap. 
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Fig. 8. — Spectra of DLA 2231—00 covering the 1656A and 1560A multiplets. Notation as in Figure 1. The best fit requires two CI 
velocity components: component 1 at V = —77 km s — 1 , or z a (, a = 2.06534, with no measurable fine structure lines, and component 2 at v 
= km s _1 , or z a (, s = 2.066122, with detected CI fine structure lines (however, N(CI**) is technically an upper limit). 
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Fig. 9. — Spectra of DLA 2231—00 low ions, v = km s 1 is located on component 2 at z a f, s = 2.066122 and marked by a vertical 
dashed blue line. Component 1 is located at v~— 77 km s — 1 . 
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Fig. 10.— J2340 CI velocity profiles over five multiplets, a) 1656A, b) 1560A, c) 1328A, d) 1280A, e) 1277A. Notation is as in Figure 
1. The components are labeled in (a), as 1 — 9 from lowest to highest relative redshift (relative to the arbitrarily chosen v = km s — 1 at 
z abs = 2.054151) and located at the following velocities: v ~ km s _1 , 13 km s — 1 , 37 km s — 1 , 44 km s — 1 , 52 km s" 1 , 55 km s -1 , 57 
km s — 1 , 83 km s — 1 , and 96 km s _1 . 
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Fig. 11. — Excitation of the fine structure level CI* (top) and CI** (bottom) caused by increasing the strength of the radiation field. We 



plot the ratio 



and 



----- ^tqy\ anc ^ nave included spontaneous radiative decay, excitation by the CMB at z a (, s =2, and collisions with 

neutral hydrogen at a temperature of T=100K. HM is the value of the Haardt-Madau background at z = 2 and is therefore a minimum 

ergs cm" 



total radiation field. Note that the radiation field must be ~> 1 X 10 18 5 ergs cm 2 s 1 Hz 1 sr 1 in order to cause significant effects 



to the level populations of the CI fine structure states at low densities (i.e. when collisions are not the dominant mechanism). 
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Fig. 12.— fl verus £2 for component 3 of DLA 0812+32, where fl = n(CI*)/n(CI) tot , f2 = n(C I**)/n(C I) tot , and n(CI) tot = n(C I) + 
n(CI*) + n(CI**). The data point is marked by a red asterisk and the la error polygon is marked in red, while the 2 a error polygon is 
blue. Theoretical tracks are indicated by black points and run from T = 10 - 10,000 K and n(H I) = 10~ 3 5 — 10 41 cm -3 . For clarity we 
have highlighted the tracks corresponding to T = 10, 100 and 1000 K in cyan. 
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Fig. 13. — n(H I) versus temperature for the allowed solutions of DLA 0812+32 component 3 from the C I analysis, la results are in 
red while 2a is in blue. 
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Fig. 14. — n(H I) versus ^^//^ for the allowed solutions of DLA 0812+32 component 3. The measured is indicated by the red 

dashed line with green dashed lines indicating the errors of ±0.2 dex. la results are in red while 2a is in blue. 
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Fig. 15. — n(H I) versus temperature for the allowed solutions of DLA 0812+32 component 3. The final solutions, as constrained by 
^jtj , are shown in black (la) and yellow (2a). 
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Fig. 16. — Example of a case, shown here for DLA 1331+17, component 1, where it is possible to constrain the radiation field using only 
the CI fine structure states and the assumption of ionization equilibrium. Log(^jy-) versus n(H I) for a selection of J v values. Solutions 
falling within the allowed range of log( ) are marked in red. It is seen that for J v total >~lxl0 -19 ergs cm -2 s _1 Hz -1 sr — 1 , there 
are no acceptable solutions. Note: This is for the Iff CI solutions and the value of the Haardt-Madau background for this object is 2.53 
xl0 — 20 ergs cm -2 s — 1 Hz -1 sr _1 . It is also apparent that the range of allowed densities is constrained to be 10 < n < 50 cm -3 . 
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Fig. 17. — Examples of the H2 absorption in DLA 0812+32. We show the relatively unblended sections here (blends are denoted by 
dotted lines) over several rotational J levels. Blue vertical dashed lines mark the three H2 components, from left to right, components 1, 
2, 3. The close association with the three CI components is clear. 
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Fig. 18. — H2 excitation diagrams for the components of DLA 0812+32. 
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Fig. 19. — Critical densities as a function of temperature for the H2 rotational J states. The subscripts u and 1 denote upper and lower 
respectively. 
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Fig. 20. — DLA 2340—00 low ions, Fe III and Ar I. Blends are indicated by orange dashed lines. The velocity space defining the three 
super-components a, b, and c are separated by vertical blue dashed lines. Relative to v = at z a i, a = 2.054151, super-component (a) is 
from —30 - 15 km s _1 , super-component (b) from 15 - 70 km s — 1 , and super-component (c) from 70 - 120 km s — 1 . For reference, the 9 
CI velocity components are indicated by red vertical tick marks. 




Fig. 21. — Examples of the H2 absorption in DLA 2340—00. We show the relatively unblended sections here (blends are denoted by 
dotted lines) over several rotational J levels. Blue vertical dashed lines mark the six H2 components we call components 1, 2, 4, 6, 8 and 9. 
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Fig. 22. — The H2 excitation diagrams for the 6 components of DLA 2340 —00 that contain H2. 
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Fig. 23. — Pressure curves and star formation rate solutions for DLA 1331+17. The blue model curves for different star formation rates 
represent logT,sFR = -co, —4, —3.5, 0.0. The CNM P m i n; P e q and Pmax solutions are indicated by the red dashed lines. Vertical, 
black dotted lines illustrate the location of the CNM stable points associated with, from left to right, P m in, Pgeo and Pmax for the case 
of background radiation only. 
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Fig. 24. — DLA 1331+17 2a range of solutions for range of CII* technique solutions. The minimal depiction CII* solution is indicated 
by the asterik, with P m i n and Pmax denoted by points at the ends of attached lines, while the maximal depletion model is denoted by a 
diamond, and the Pmin and Pmax at the ends of dashed lines. For comparison, the full range of CI solutions arc given as a grouping of 
points in the upper right. It is clear that these fall at higher pressures and densities than than the CII* results. 
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Fig. 25. — Summary of all DLAs, compared with the results of the C II* technique (diamonds). The 2a CI results are the shaded regions, 
and the best-fit (minimized x 2 ) for each is marked by a filled circle. It is clear that the CI results are consistently higher in density and 
therefore pressure. Note, some regions/solutions are overlapping. The H2 derived pressures are denoted by light grey crosses and asterisks 
for components in DLA 0812+32 and DLA 2340—00 respectively. The solid black line denotes T = 10 K. The resulting median values for 



the C I sample are: <n(H I)> 
a T = 52 K, and <r log(P/k) 
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139 K, and <log(P/k)> = 2.57 cm 3 K , with standard deviations a. 
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Fig. 26. — The resultant cloud size in pc (top axis) as a function of the fraction of CII, f(CII), associated with the CI cloud. Here, f(CII 
) = 1 means that 100% of the measured N(CII) is associated with the CI cloud. In this case, for DLA 0812+32, component 3, log( c ^ ) 
= 3.1 ± 0.2. As the fraction of associated CII decreases, the n(H I) increases and hence, the resultant cloud size decreases. 
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TABLE 1 
DETAILS OF OBSERVATIONS 



Quasar 


Telescope/Instrument 


Date 


Resolution 


Total Exp. Time 








FWHM [km s" 1 ] 


N 


FJ0812+32 


Keck/HIRES 


Mar. 16/17, 2005 


8.33 


14400 


FJ0812+32 


Keck/HIRES 


Jan. 12, 2008 


6.25 


7200 


Q1331+17 


Keck/HIRES 


Jan. 10, 1994 


6.25 


3600 


Q1331+17 


VLT/UVES 


Dec. 18, 2001 


7.0 


4500 


J2100-06 


Keck/HIRES 


Sept. 18, 2007 


6.25 


10800 


Q2231-00 


Keck/HIRES 


Jan. 11, 1995 


8.33 


5400 


J2340-00 


Keck/HIRES 


Aug. 18/19, 2006 


6.25 


15000 



TABLE 2 
C I Data 



Quasar 


cmp z aba FWHM.^t" 
[km s — 1 ] 


b 

[km s _1 ] 


logN(CI) 
[cm -2 ] 


logN(CI*) 
[cm" 2 ] 


logN(CI**)] 
[cm" 2 ] 


rp()l b 

ex 

[K] 


(fl, f2) c 


FJ0812+32 


1 2.066780(7) 6.25/8.33 


1.20 ± 0.14 


12.96 ± 0.04 


12.92 ± 0.02 


12.35 ± 0.04 


19.8 


(0.43,0.11) 





FJ0812+32 1 2.625808(2) 6.25/8.33 3.25 ± 1.00 

FJ0812+32 2 2.6263247(8) ... 2.57 ± 0.56 

FJ0812+32 3 2.626491(1) ... 0.33 e ± 0.05 



12.13 ± 0.05 11.68 ± 0.16 11.37 ± 0.270 11.1 (0.23, 0.11) 
12.70 ± 0.02 12.32 ± 0.04 <12.39(1<t) 12.0 (0.22, 0.26) d 
13.30 ± 0.23 13.02 ± 0.03 12.47 ± 0.05 13.5 (0.31,0.09) 



Q1331+17 1 1.7763702(9) 6.25/7.0 5.08 ± 0.24 13.08 ± 0.02 12.59 ± 0.02 11.90 ± 0.14 10.6 (0.23, 0.05) 

Q1331+17 2 1.776524(1) ... 0.55 ± 0.13 13.06 ± 0.13 12.05 ± 0.07 <11.7 6.9 f 

Q1331+17 3 1.77664(5) ... 24.45 ± 6.2 12.51 ± 0.11 <12.2 <12.25 <13.0 



Q2100-06 1 3.089776(7) 6.25 0.20 ± 0.3 12.28 ± 0.21 12.01 ± 0.13 11.56 ± 0.28 13.7 (0.31, 0.11) 

Q2100-06 2 3.091463(4) ... 3.89 ± 0.69 12.57 ± 0.03 12.31 ± 0.08 13.9 

Q2100-06 3 3.09236(2) ... 12.77 ± 1.84 12.61 ± 0.05 12.16 ± 0.12 11.1 



Q2231-00 1 2.06534(3) 8.33 19.08 ± 4.37 12.55 ± 0.08 

Q2231-00 2 2.066122(5) ... 5.36 ± 1.17 12.65 ± 0.05 12.41 ± 0.09 11.60 ± 0.16 14.3 (0.35, 0.05) 



J2340- 


-00 


1 


2.054151(4) 


6.25 3.36 ± 0.80 


12.32 ± 0.04 


11.96 ± 


0.12 




12.2 




J2340- 


-00 


2 


2.054285(3) 


0.24 ± 0.12 


12.54 ± 0.17 


12.28 ± 


0.09 


11.84 ± 0.14 


13.9 


(0.32,0.11) 


J2340- 


-00 


3 


2.054526(1) 


2.55± 0.16 


13.33 ± 0.02 


12.98 ± 


0.02 


11.52 ± 0.33 


12.4 


(0.31, 0.01) 


J2340- 


-00 


4 


2.054602(3) 


1.39 ± 0.70 


12.97 ± 0.12 


12.53 ± 


0.12 


12.12 ± 0.13 


11.2 


(0.24, 0.09) 


J2340- 


-00 


5 


2.05469(3) 


7.01 ± 3.21 


13.03 ± 0.08 


12.46 ± 


0.31 


12.02 ± 0.34 


9.8 


(0.20, 0.07) 


J2340- 


-00 


6 


2.054711(9) 


1.19 ± 2.31 


12.54 ± 0.67 


12.64 ± 


0.36 


11.96 ± 0.26 


27.2 


(0.50, 0.11) 


J2340- 


-00 


7 


2.054727(5) 


0.61 ± 0.29 


13.48 ± 0.10 


12.59 ± 


0.37 




7.5 




J2340- 


-00 


8 


2.05499(1) 


9.12 ± 2.11 


12.56 ± 0.07 


11.96 ± 


0.16 


11.70 ± 0.29 


9.5 


(0.18, 0.10) 


J2340- 


-00 


9 


2.055131(4) 


2.92 ± 0.90 


12.47 ± 0.07 


12.15 ± 


0.08 


11.72 ± 0.16 


12.9 


(0.29, 0.11) 



J2340-00 1 1.36027 6.25 1.58 ± 12.83 ± 0.07 12.28 ± 0.07 11.96 ± 0.14 10.0 (0.20, 0.10) 

J2340-00 2 1.36049 ... 1.03 ± 12.68 ± 0.08 11.78 ± 0.29 7.4 

J2340-00 3 1.36061 ... 9.75± 11.74 ± 0.27 12.42 ± 0.15 11.96 ± 0.26 ... g (0.64, 0.22) 

J2340-00 4 1.36088 ... 0.64 ± 13.23 ± 0.28 12.77 ± 0.08 11.81 ± 0.21 10.9 (0.25, 0.03) 



a Instrumental FWHM 

Tei i s the excitation temperature derived from the CI and CI* states. This should be consistent with and generally is higher than Tcmb- 
c fl and f2 as defined in the text 

^ (fl, f2 ) derived using the upper limit on N(C I**) 

e Dopplcr parameter determined from the curve of growth method in [jorgenson ct al. (2009) 

^ T GX < Tcmb- This unphysical result is likely due to measurement error. Indeed, accounting for the given error, T eK — 8 K consistent with Tcmb 
= 7.6 K. 

g N(CI*) >> N(CI) which implies T c:c = -50 K. 



TABLE 3 

Other Relevant DLA Data 



Quasar 


■^abs 


logN(H I) 

[cm- 2 ] 


logN(cir) 

[cm -2 ] 


ic 

[ergs s~ ^-H - *] 


[M/H] a 


[Fe/H] 6 


logN(CT ot ) 


log (%¥) c 

[cm- 2 ] 


Av 

[km s- 1 ] 


FJ0812+32 


2.066780 


21.50 ± 0.20 


13.62 d 


-27.4 d 


-1.83 ± 0.20 


-2.06 ± 0.02 


13.30 


2.57 ± 0.20 


26 








14.42 e 


-26.6 e 












FJ0812+32 


2.62633 


21.35 ± 0.10 


14.30 ± 0.01 


-26.56 ± 0.10 


-O.Sl f ± 0.10 


-1.629 ± 0.10 


13.63 


3.10 ± 0.20 


70 


Q1331+17 


1.77636 


21.14 ± 0.08 


<13.56' 1 


<-27.16 


-1.37 ± 0.01 


-1.97 ± 0.01 


13.56 


2.40 ± 0.20 


72 


Q2100-06 


3.09237 


21.05 ± 0.15 


14.09 ± 0.01 


-26.48 ± 0.15 


-0.73 ±0.15 


-1.20 ± 0.02 


13.17 


3.39 ± 0.02 


187 


Q2231-00 


2.06615 


20.56 ± 0.10 


13.71 ± 0.04 


-26.37 ± 0.11 


-0.74 ±0.16 


-1.40 ±0.07 


13.04 


2.88 ± 0.20 


122 


J2340-00 


2.05452 


20.35 ± 0.15 


13.84 ± 0.04 


-26.03 ± 0.15 


-0.74 ±0.16 


-0.92 ± 0.03 


14.09 


1.86 ± 0.20 


104 



a Logarithmic ct-mctal abundance with respect to solar, determined using Si II unless otherwise labeled 

Logarithmic Fc abundance with respect to solar 
c N(CII) calculated assuming the "minimal depletion model" from I Wolfe et al. I J2004D where [C/H] = [Si/H] - 0.2. Sec text for details. 

Assume low cool DLA with average low cool i c . (The measured upper limit on log N(CII*) — 14.06 cm"" 2 , which includes blending) 
c Assume high cool DLA with average high cool l c 
^ Global metallicity based on Zn II measurement. 
s Based on Cr II instead of Fe II 

k N(CII*) as determined by fitting constrained by Si II components. Blending means that this is still technically an upper limit. 



Jorgenson et al. 



TABLE 4 
CII* TECHNIQUE SOLUTIONS 



Quasar 


z abs 


local a 

Ji> ° Ca a 


7 Bkdd 
■'v 






/lxlO- 19 


/lxlO- 20 


r J Do lZ-\-oZ 


z.UDo ( o 


U.4o 


Z.OO 


FJ0812+32 


2.06678 


3.46 c 


2.56 


FJ0812+32 


2.06678 


36.38 d 


2.56 


FJ0812+329 io!,ai 


2.62633 


7.17 


2.45 


FJ0812+32^ mpl 


2.62633 


19.2 


2.45 


FJ0812+32^ p2 


2.62633 


6.3 


2.45 


FJ0812+32^ p3 


2.62633 


259.9 


2.45 


Q1331+17 


1.77636 


3.09 


2.53 


J2100-06 


3.09237 


17 


2.33 


Q2231-00 


2.06615 


24.7 


2.56 


J2340-00 


2.05452 


52.4 e ,10.8 / 


2.56 



a Evaluated at 1500A; Error is approximately ± 50% 

b CII* in forest. Assumed low £ c =-27 A and [Fe/Si]=0.0, 

the maximal depletion model 

c C II* in forest. Assumed low £ c =-27 A and [Fe/Si]=-0.2, 
the minimal depletion model 

d C II* in forest. Assumed high £ c =-26.6 and [Fe/Si] = -0.2 
° J v local determined for components 3, 4, 5 
f j^iocai determined for components 6 & 7 



TABLE 5 
CI TECHNIQUE SOLUTIONS 



DLA, comp. 


j total 


l°g(%T-) 




1(7 (2(7) C I constraints 








best x 2 


fit 


I a 




x'io- 19 


n(H I ) 




T 




log(P/k) 


n(H I ) 


T 


log(P/k) 










[cm' 3 ] 




[K] 




[cm~ 3 K] 


[cm- 3 ] 


[K] 


[cm~ 3 K] 


[PC] 


0812+32^=2.06, 1(a) 


3.72* 


2.57+0.2 


302 - 1096(302 - 1202) 


32 


— JU^Oi — (JO ) 


4.18 


- 4.54(4.10 - 4.58) 


316 


50 


4.20 


3.2 


0812+32z o[)s =2.06, 1(b) 


0.69 c 


2.77+0.2 


40 - 83(35 - 91) 


251 — 


1000(200 — 1259) 


4.26 


- 4.66(4.08 - 4.80) 


79 


251 


4.30 


13.0 


0812+32^=2. 06 , 1(c) 


36.64 d 


2.57+0.2 


- (3802 - 12589) 




— ^zu — ZO J 




- (4.88 - 5.40) 


12589 


20 


5.40 


0.1 


0812+32 iaiis= 2.62 9loi ' ai ,l 


7.4 


3.10+0.2 


72 - 120(23 - 417) 


32 


- 50(10 - 126) 


3.48 


- 3.58(2.36 - 4.34) 


76 


40 


3.49 


9.5 




11. 0" 


3.10+0.2 


105 - 151(46 - 602) 


25 


- 32(10 - 79) 


3.52 


- 3.58(2.66 - 4.38) 












3.8^ 


3.10+0.2 


50 — 138(13 — 275) 


32 


— 53(10 — 200) 


3.50 


— 3.64(2.12 — 4.42) 










0812+32 iQiis= 2.i32 9loi ' a ',3 


7.4 


3.10+0.2 


87 - 398(72 - 549) 


32 - 


- 126(25 — 251) 


3.74 


- 4.30(3.54 - 4.80) 


100 


79 


3.90 


7.2 


11.0' 


3.10+0.2 


138 - 602(105 - 912) 


25 


- 79(20 - 158) 


3.74 


- 4.38(3.52 - 4.76) 












3.8 / 


3.10+0.2 


66 - 275(55 - 417) 


40 - 


- 158(32 - 398) 


3.74 


- 4.28(3.54 - 4.90) 










0812+32 Zab3= 2.62 dt, ',l 


19.4 


3.71+0.2 


14 — 57(5 — 158) 


25 


- 79(10 - 126) 


2.86 


— 3.32(1.66 — 4.60) 


14 


50 


2.86 


5.9 




29. e 


3.71+0.2 


8 - 66(8 - 209) 


10 


- 32(10 - 251) 


1.92 


- 3.12(1.92 - 4.40) 












9.9 / 


3.71+0.2 


14 - 42(2 - 105) 


63 — 




3.32 


- 3.80(1.22 - 4.86) 










0812+32,^^=2. 62 d "\3 


260 


3.06+0.2 


1513-12020(1513-12587) 


10 


- 20(10 - 32) 


4.18 


- 5.38(4.18 - 5.60) 


3161 


10 


4.50 


0.1 




390.1 e 


3.06+0.2 


2629-12587(2187-12587) 


10 


- 20(10 - 32) 


4.42 


- 5.40(4.34 - 5.60) 












130 2^ 


o nfi+n 9 

O.UOxU.Z 




10 


- 20(10 - 40) 


3 88 


K 88 ^ 49*1 
— O.UO^O.OO — 










1331 + 17,1 


3.3 


2.40±0.2 


- (229 - 457) 




- (16 - 20) 




- (3.66 - 3.86) 


331 


16 


3.72 


1.4 




4.9" 


2.40±0.2 


- (417 - 1738) 




- (13 - 16) 




- (3.82 - 4.34) 












1.8 / 


9 4.n_|_n 9 






- (16 - 32) 




— 4n — ^ 










ziuu — UD, 1 


17 2 


^ ^Q-l-n on 
o . o y It u . z u 


99 ^9^/99 1 9^S7A 


10 - 


251(10 - 10000) 


2 34 


A 1.A in\ 

— 4.(0(Z.04 — O.IUJ 


fid 
oy 


ou 


^ X.A 


0.0 




25 7 e 


^ ^Q-l~n on 


AA QVlCAA 1 9^S7A 

44 — Of — IZOOf ) 


10 - 


158(10 - 10000) 


2 64 


A 7 A (O P.A Q id) 
— 4. |4^Z.D4 — O.IUJ 












8.7 


3.39±0.20 


24 — 347(10 — 7584) 


20 - 


501(10 - 10000) 


2.88 


— 4.94(2.02 — 7.88J 










2231-00,2 


25.0 


2.88±0.2 


199 - 1513(174 - 1819) 


13 


- 25(10 - 32) 


3.40 


- 4.48(3.24 - 4.66) 


363 


20 


3.86 


0.3 






Z.OO±U.Z 


o4 1 — z0oU(o4f — oUiyj 


10 


- 20(10 - 25) 


0.04 


— 4.0z(0.D4 — 4.(0) 












12.6^ 


2.88+0.2 


151 — 758(79 — 832) 


16 


— 40(10 — 50) 


3.54 


— 4.28(2.90 — 4.42) 










2340 — 00,2 


0.39 s 


1.86+0.2 


48 — 240(30 — 525) 


63 - 


501(40 - 1585) 


3.84 


— 4.80(3.38 — 5.68) 


66 


200 


4.12 


1.1 


2340 — 00,3 
























Zo4U — U U,4 






/is i on/ 1 ^ ^ &\ 

4o — ±ZU(00 — OlDJ 


63 - 


- 158(40 - 794) 


3 78 


— o.yo(o.4z — O.IZJ 


A 8 


1 ^8 


O.OO 


1.0 


934n nn ^ 
zoiu — uu,o 


„ 


„ 


99 1 90(7 41 7\ 


25 - 


158(10 - 1259) 


3 04 


4 nnfl 89 ^ 49"! 
— 4.UUI l.OZ — G.4ZJ 


35 


100 


3 54 


2 1 


9^4n nn s. 


„ 




Q1 ('7 1 F,fi\ 
oo — y±^f — iod j 


40 - 


- 100(10 - 316) 


3 44 


Q ftO?1 S9 4 4n^ 

— O.DZ^l.OZ — 4.4U ) 


35 


79 


3 44 


2 1 


9^4(1 — no 9 


„ 


„ 


40 _ 1 ^8^29 — 275^1 


63 - 


- 200(32 - 631) 


3.82 


— 4 1 Of 3 30 — 4 98^ 


48 


158 


3.88 


1.5 


2340-00,2 


_h 


2.24+0.2 




















2340-00,3 


52.7 


1.74+0.2 




















2340-00,4 


52.7 


1.74+0.2 




















2340-00,5 


52.7 


1.74+0.2 




















2340-00,8 


11.0 


2.34+0.2 


- (661 - 5011) 




- (10 - 20) 




- (3.82 - 5.0) 


661 


10 


3.82 


0.03 




14. 5 C 


2.34+0.2 


- (1318 - 7943) 




- (10 - 20) 




- (4.12 - 5.10) 












5.7 / 


2.34+0.2 


- (219 - 1819) 




- (10 - 32) 




- (3.34 - 4.56) 










2340-00,9 


11 


2.34+0.2 


- (912 - 6606) 




- (13 - 40) 




- (4.12 - 5.32) 


912 


16 


4.16 


0.02 




14. 5 e 


2.34+0.2 


- (1513 - 12588) 




- (13 - 40) 




- (4.28 - 5.60) 












5.7 / 


2.34+0.2 


- (363 - 2884) 




- (16 - 63) 




- (3.80 - 5.06) 











a This estimated cloud size, £ — N(H I )/n(H I ), is technically an upper limit. 

The most likely model: Minimal depletion model, low i c object with [Fc/Si] — —0.2 
c Maximal depletion model, assuming the low cool £ c with [Fe/Si] — 0.0 
d Minimal depletion model, assuming the high cool i c with [Fe/Si] — —0.2 
c Includes +50% error on jJ ocal 
f Includes -50% error on J u local 

s Minimum background assumption J v ~ Haardt-Madau for this object 
h Solar (or super-solar) Fe/H precludes calculation. 



TABLE 6 

C I Technique solutions: full J v total grid 



DLA 


corrip. 


CI > 


n(H I)(2ct) 


T(2a) 


logfP Mi 2a) 


Jvtotal /1x10 _19 (2<t) 


Consistent? ^ 






[cm -3 ] 


[K] 


[cm- 3 K] 


a 




0812+32 Zabs=2 .06 


i 


2.57±0.2 


>36 (>32) 


20-1258 (13-1585) 


4.08-5.42 (4.0-5.50) 


0.41-195 (0.41-275) 


yes 


0812+32 Zabs=2 . 62 f io6ai 


i 


3.10±0.2 


>6 (>0.1) 


<4000(-) 


3.04-5.10 (0.08- 5.70) 


<773 (<773) 


yes 


U81^+cS/,z ol)s= 2.62 M 


o 
o 


o i n-i-n 9 

O. ±UlCU.Z 




SDoUU^- ) 


'i ^ ac\ fi in ^ ^^^^ 


_^ ( / O < > OJ 


yes 


0812+32 Zabs= 2.62 d ™ 


1 


3.71±0.2 


2-57(0.002-4166) 


- (-) 


1.72-4.68 (<5.34) 


<34 (<1082) 


yes 


0812+32 z . - 2 6 2 div 

^abs — ■i-u^i 


3 


3.06±0.2 


>7 (>6) 


<7943(-) 


3.40-5.40 (3.18- 5.60) 


<923 (<923) 


yes 


1331+17 


1 


2.40 ±0.2 


11-44 (10-12021) 


79-794 (<794) 


3.50-4.04 (3.20-5.08) 


0.34-0.94 (0.34-86) 


yes 


2100—06 


1 


3.39±0.20 


>2 (>0.03) 




2.14—5.60 (-) 






2231-00 


2 


2.88±0.2 


>6 (>3) 


<6300 (<7950) 


3.40-5.30 (2.54-5.40) 


<751 (<751) 


yes 


2340-00 


2 


1.86±0.2 


>48 (>30) 


20-500 (16-1585) 


3.82-5.50 (3.38-5.90) 


<19 (<21) 


_ 


2340-00 


3 




>132 (>110) 


13-32 (13-32) 


3.62-5.20 (3.54-5.20) 


0.92-67 (0.44-24) 


no 


2340-00 


4 


" 


48-190 (>33) 


40-158 (16-794) 


3.78-3.98 (3.42-5.68) 


0.39-1.1 (0.39-21) 


no 


2340-00 


5 




>22 (>7) 


13-158 (<1259) 


3.04-5.20 (1.82-5.80) 


0.39-23.8 (0.39-26.7) 


no 


2340-00 


8 




35-120 (>7) 


32-100 (<316) 


3.44-3.62 (1.82-5.30) 


0.39-0.92 (0.39-26.7) 


yes 


2340-00 


9 




47-316 (>29) 


32-200 (13-631) 


3.82-4.10 (3.30-5.60) 


0.39-1.4 (0.39-23.8) 


yes 


2340-00 


2 


2.24±0.2 


>24 (>13) 


20-1995 (13-6310) 


3.82-5.50 (3.38-5.90) 


0.39-53(0.39-66.7) 




2340-00 


3 


1.74±0.2 


>132 (>110) 


13-32 (13-32) 


3.62-5.20 (3.54-5.20) 


0.92-67 (0.39-24) 




2340-00 


4 


1.74±0.2 


60-190 (>44) 


40-126 (16-501) 


3.78-3.94 (3.44-5.68) 


<1.1 (<21) 




2340-00 


5 


1.74±0.2 


>32 (>12) 


13-126 (<794) 


3.10-5.20 (2.06-5.80) 


<23.8 (<26.7) 




2340-00 


8 


2.34±0.2 


34-120 (>7) 


32-100 (<316) 


3.44-3.62 (1.82-5.30) 


<9.2 (<26.7) 




2340-00 


9 


2.34±0.2 


48-316 (>29) 


32-200 (13-631) 


3.82-4.10 (3.30-5.60) 


<1.4 (<23.8) 





a [ergs cm 2 s 1 Hz 1 sr 1 ] 

b Is the radiation field derived via the C II* technique consistent with the range allowed by the C I data? 



Using C I fine structure to probe physical conditions in DLAs 



TABLE 7 

Dust to Gas Ratio Component Analysis of DLA 0812+32 



comp. 1 


comp. 2 


comp. 3 


comp. 4 


z abs =2. 625890 


z abs =2.626310 


z ai)a =2.626491 


z a6s =2.626447 


b = 18.99±1.91 [km s" 1 ] 


b = 6.23±0.35 [km s" 1 ] 


b = 0.33±0.05 a [km s" 1 ] 


b = 5.6±3.3 [km s" 1 ] 



1 o Qi _|_ n n/i 
Iz.yl it U.U4 


iq ic J. n no 
lo.lO it U.UZ 


\ 1 1. / 1 1CT J 


i o c;o -U n D7 
iz.oz it U.U f 


i o A a _i_ n n a 
VAAt± ± U.U4 


12.yo ± U.U2 


13.00 & 


TO QO _L C\ O A 
1Z.OZ it U.Z4 


13.00 ± 0.2 


13.69 ± 0.01 


15.14 ± 0.17 


13.08 ± 0.05 


20.41 


20.93 


20.97 


20.29 


—U.oo ± U.U4 


—U.oo ± U.Uz 


—0.58 


n e; o _i_ n o a 
— U.oo ± U.24 


i i q _i_ n 04. 
— ±.J.O m u.u^ 


— X.^tJL m U.UZ 


— z.oy 


1 40 4- 07 

— 1 .^U HI U.U I 


_n crc _i_ n 
— u.uu m u.uu 


_n oq _|_ n n*3 
— u.oo m u.uo 


<— 2.31 


_n QO -j- fl 19 

U.O.^ HI U.-L^ 


0.09 


0.13 


~0.16 


0.13 


-1.04 


-0.90 


-0.78 


-0.90 


3.71 ± 0.2 


3.56 ± 0.2 


3.06 ± 0.2 




-26.37 


-26.76 


-25.35 


-26.73 


19.2 


6.3 


259.9 


2.5 


2.34xl0" 5 


9.0xl0~ 3 


0.14 




3.53X10- 11 


6.74X10" 14 


1.53X10" 14 




5.97X10- 11 


8.96xl0~ 14 






0.08 


3.6xl0~ 4 


3.9X10" 5 




0.98 


0.91 


0.88 




0.08 


3.3xl0~ 4 


3.41xl0~ 5 




102 


64 


47 




0.37 


0.17 


0.36 




21 


11 


37 





N(Cr II) [cm" 2 ] 
N(Zn II) [cm" 2 ] 
N(CII*) [cm" 2 ] 
N(H I) d [cm" 2 ] 

[Zn/H] 

[Cr/H] 

[Cr/Zn] 

K e 

logi K 
CII 

T c 

Jjocal /IQ-Wf 
f H 2 

fl) [s- 1 ] 
fl [s- 1 ] 

[K] 

j^/io- 19c 

n(H 1)9 [cm" 3 ] 



a Dopplcr parameter fixed to match that of C I as determined by IJorgenson et all (2009) 

N(Zn II) fixed by upper limit on N(0 I) assuming solar relative abundances. See text for details. 
c Summed over the two components required by VPFIT, as explained in the text. 
d N(H I) scaled to trace N(Zn II) 

c Dust to gas ratio relative to Milky Way, as defined in the text. 

[ergs cm" 2 s _1 Hz -1 sr _1 ] 
s Density derived from the H2 as explained in the text. 



Jorgenson et al. 



TABLE 8 

DLA 0812+32 Molecular Hydrogen 



comp 


ion 


z 


z abs 


b 

[km s" 1 ] 




log N 
[cm- 2 l 

L J 


(Tins, AT 
i og jv 


T _j a 

L J 


comp 1 




















H2J0 


2.625812 


0.000011 


1.03 


0.11 


15.03 


0.15 


- 




H2J1 


" 


" 


" 


" 


15.26 


0.15 


102 ±\ 




H2J2 


" 


" 


" 


" 


14.04 


0.06 


132 +1 




H2J3 


" 


" 


" 


" 


13.39 


0.10 


150 + _ 2 2 




H2J4 


" 


" 


" 


" 


<12.64 


- 


222 +1° 




H2J5 


" 


" 


" 


" 


12.22 


0.61 


257 tH 




total N(H 2 ) = 


3.01 xlO 15 














log (total N(H 2 )) = 


15.48 
















f b > 


2.69xl0~ 6 














comp 2 




















H2J0 


2.626326 


0.000001 


2.67 


0.09 


<18.45 c 








H2J1 










18.25 


0.04 


~ . -4-9 A 

64 tit 




H2J2 










15.49 


0.06 


61 t\ 




H2J3 










14.43 


0.01 


83 t\ 




H2J4 


" 


" 




" 


13.18 


0.10 


119 t\ 




H2J5 


" 




" 


" 


<12.60 


- 


151 t\ 




total N(H 2 ) < 


4.61X10 18 














log (total N(H 2 )) < 


18.66 
















{ b > 


4.1 xl0~ 3 














comp 3, model l d 




















H2J0 


2.626491 


0.000001 


0.81 d 


- 


19.79 


0.03 


- 




H2J1 


)5 




55 




19.15 


0.03 


47 +\ 




H2J2 


55 


55 


55 




16.60 


0.03 


57 +\ 




H2J3 


55 


55 


55 




15.11 


0.05 


74 t\ 




H2J4 


55 


55 


55 




13.99 


0.04 


110 




total N(H 2 ) = 


7.58X10 19 
















log (total N(H 2 )) = 


19.88 
















f b > 


0.06 














comp 3, model 2 e 




















H2J0 


2.626494 


0.000001 


1.19 e 




19.81 


0.03 






H2J1 


)5 




55 




19.13 


0.03 


45 +\ 




H2J2 


55 


55 


55 




16.21 


0.11 


52 + 2 




H2J3 


55 


55 


55 




14.88 


0.07 


71 t\ 




H2J4 


55 


55 


55 




13.97 


0.04 


109 t J 



total N(H 2 ) = 
log (total N(H 2 )) 



7.81 xlO 19 
19.89 
0.07 



a Excitation temperature between rotational level J and J — 0. 

k Molecular fraction f calculated using N(H l^ total — 21.35 cm -2 , so that these values are technically lower limits. 

c We report the 2cr upper limit because blending with the stronger J— line of component 3 make the formal errors large. 

^ Model 1: Doppler parameter tied to that of C I. Note that the J— and J— 1 transitions are heavily saturated and therefore 

the resultant log N is determined by the damping wings and is insensitive to the choice of b. 

c Model 2: Doppler parameter determined by the J— 3 state. 



Using C I fine structure to probe physical conditions in DLAs 



TABLE 9 

AODM Component Analysis of DLA 2340-00 



Super-Component a b c 

Av 90 -30-15 km s- 1 15-70 km s" 1 70-120 km s" 1 



C I Component 






(8,9) 


N(Fc II) [cm" 2 ] 


14.55 ± 0.05 


14.63 ± 0.04 


14.06 ± 0.01 


N(S II) [cm" 2 ] 


14.17 ± 0.01 


14.73 ± 0.01 


14.31 ± 0.01 


N(Cir) [cm" 2 ] 


12.53 ± 0.09 


13.60 ± 0.01 


12.99 ± 0.04 


N(H I) a [cm- 2 ] 


19.58 


20.13 


19.72 


[b e/H\ 


n a 7 _l_ n 
— U.4f ± U.Uo 


— U.yo ± U.U4 


— 1.1U ± U.U1 


[S/H] 


-0.56 ± 0.01 


-0.56 ± 0.01 


-0.56 ± 0.01 


[Fe/Met] 


0.09 ± 0.05 


-0.39 ± 0.04 


-0.54 ± 0.01 




-0.06 


0.06 


0.09 


logloK 


_ 


-1.21 


-1.03 


N(CI)* otai 


12.96 


14.01 


13.00 


lc 


2.24 


1.74 


2.34 


-26.56 


-26.05 


-26.24 


Jjocal /1Q-19C 


_d 


52.4 


10.8 


N(Fe III) [cm- 2 ] 


<13.87 e 


< 13.94 e 


<13.89 e 


Fe III/Fc Il-f 


<0.21 


<0.21 


<0.66 


N(Ar I) [cm- 2 ] 


13.31 


13.83 


13.65 


[Ar/S] 9 


0.12 


0.08 


0.32 


N(Ni II) [cm- 2 ] 


13.21 ± 0.02 


13.50 ± 0.01 


13.14 ± 0.02 


[Ni/H] 


-0.55 ± 0.02 


-0.82 ± 0.01 


-0.79 ± 0.02 


[N/Met] 


0.01 ± 0.02 


- 0.26 ± 0.01 


- 0.23 ± 0.03 


n 


-0.01 


0.02 


0.01 


logK 




-1.67 


- 1.97 



a N(H I) scaled to trace N(S II) 

b Dust to gas ratio relative to Milky Way, defined in the text. Here we have used 

S II instead of Si II. 

c [ergs cm -2 s -1 Hz -1 sr -1 ] 

d We did not determine J L , local for this component because the super-solar Fe II 
measurement. 

c AODM measurements taken as upper limits because of possible blending with 
the forest. 

f Fe III/Fc II ~0.3 means partially ionized, HI/H = 0.5 

S If [Ar/S] > -0.2 then x <0.1 (but low Ar/S docs not require x>>0). 
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DLA 2340 



TABLE 10 
-00 Molecular Hydrogen 



component™ 


ion 


z abs 




b 

[km s- 1 ] 




log N 
[cm" 2 ] 


VlogN 


T -j 6 
[K] 


comp 1 


H2J0 


2.054165 


0.000001 


2.31 


0.07 


15.262 


0.040 






H2J1 


2.054165 




2.31 


0.00 


15.94 


0.049 


266 ±f 




H2J2 


2.054165 




2.31 


0.00 


14.904 


0.065 


211 +\ 




H2J3 


2.054165 




2.31 


0.00 


14.251 


0.049 


191 ±\ 




H2J4 


2.054165 




2.31 


0.00 


<13.39 


- 


262 t\ 




H2J5 


2.054165 




2.31 


0.00 


<12.90 


— 


287 ±1 




log (total N(H 2 )) = 
f c = 


16.06 
1.03xl0~ 4 














comp 2 


H2J0 


2.054291 


0.000001 


1.25 


0.09 


14.657 


0.073 






H2J1 


2.054291 




1.25 


0.00 


15.292 


0.079 


232 ±\ 




H2J2 


2.054291 




1.25 


0.00 


14.478 


0.079 


253 ±1 




H2J3 


2.054291 




1.25 


0.00 


14.299 


0.059 


265 ±1 




H2J4 


2.054291 




1.25 


0.00 


<13.08 


- 


293 + l 




H2J5 


2.054291 




1.25 


0.00 


<13.23 


_ 


378 ±1° 




log (total N(H 2 )) = 
f c = 


15.47 
2.63 xl0~ 5 














comp 4 


H2J0 


2.054573 


0.000002 


4.62 


0.16 


17.269 


0.082 






H2J1 


2.054573 




4.62 


0.00 


17.955 


0.054 


276 + 2? 

— ZD 




H2J2 


2.054573 




4.62 


0.00 


17.045 


0.129 


241 +]l 
— 12 




H2J3 


2.054573 




4.62 


0.00 


15.326 


0.028 


136 1?, 




H2J4 


2.054573 




4.62 


0.00 


13.812 


0.041 


168 l 2 . 




H2J5 


2.054573 




4.62 


0.00 


13.317 


0.139 


203 t 2 




log (total N(H 2 )) = 
f c = 


18.08 
1. 06X10" 2 














comp 6 


H2J0 


2.054714 


0.000001 


5.06 


0.11 


15.997 


0.047 






H2J1 


2.054714 




5.06 


0.00 


16.825 


0.065 


587 tf 4 




H2J2 


2.054714 




5.06 


0.00 


16.059 


0.039 


349 ±1 




H2J3 


2.054714 




5.06 


0.00 


15.836 


0.019 


300 ±1 




H2J4 


2.054714 




5.06 


0.00 


14.355 


0.015 


286 ±1 




H2J5 


2.054714 




5.06 


0.00 


13.932 


0.040 


310 ±\ 




log (total N(H 2 ) = 
f c = 


16.98 
8.51xl0~ 4 














comp 8 


H2J0 


2.054986 


0.000001 


3.64 


0.10 


15.76 


0.042 






H2J1 


2.054986 




3.64 


0.00 


16.55 


0.067 


475 tf 7 




H2J2 


2.054986 




3.64 


0.00 


15.63 


0.047 


271 +J 2 




H2J3 


2.054986 




3.64 


0.00 


15.13 


0.027 


229 +2 




H2J4 


2.054986 




3.64 


0.00 


13.75 


0.046 


251 +\ 




H2J5 


2.054986 




3.64 


0.00 


11.90 


3.520 


207 ±|f° 




log (total N(H 2 ) = 
f c = 


16.67 
4.21X10" 4 












comp 9 


H2J0 


2.055135 


0.000001 


1.80 


0.07 


16.735 


0.074 






H2J1 


2.055135 




1.80 


0.00 


17.200 


0.057 


151 tl 




H2J2 


2.055135 




1.80 


0.00 


16.035 


0.088 


159 +2 




H2J3 


2.055135 




1.80 


0.00 


14.753 


0.049 


135 ±\ 




H2J4 


2.055135 




1.80 


0.00 


<13.10 




162 +1 




H2J5 


2.055135 




1.80 


0.00 


<13.02 




212 +1 




log (total N(H 2 ) = 
f c = 


17.35 
2.00xl0~ 3 















a Excitation temperature between rotational level J and J — 0. 
k Components are numbered by the closest z C I component 
c f calculated assuming logN(HI) — 20.35 
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TABLE 11 

H 2 Component Analysis of DLA 2340-00 - using AODM values 



Super-Component 


a 




a 




b 




b 




c 




c 




component 


l a 




2 a 




4 




6 




8 




9 






6.04x10' 


-4 


1.55x10" 


-4 


1.74x10" 


-2 


1.41x10- 


-3 


1.79x10- 


-3 


8.47x10" 


-3 


[K] 


266 




232 




276 




587 




475 




151 




A> Is" 1 ] 


9.44x10" 


-11 


1.84x10" 


■10 


2.43x10" 


12 


1.36x10" 


10 


6.31x10" 


11 


1.84x10" 


12 


01 [s- 1 ] 


1.00x10" 


■10 


8.42x10" 


■10 


2.52x10" 


12 


1.44x10" 


10 


2.62x10" 


12 


5.89x10" 


12 


Sself 


2.85x10" 


-2 


7.92x10" 


-2 


8.72x10" 


-4 


5.83x10- 


-3 


9.88x10- 


-3 


3.07x10" 


-3 


^dust 










0.993 




0.993 




0.996 




0.996 




Stotal 










8.66x10- 


-4 


5.78x10- 


-3 


9.84x10- 


-3 


3.06x10" 


-3 


j^ w /io- 19c 










2.25 




18.84 




5.13 




< 0.48 




n(H I) d [cm" 3 ] 










1629 




10509 




3595 




377 





a Super-solar Fe II measurement, therefore k was not sensible 

b f is calculated using the N(H I) of the Super-Component, i.e. N(H I) = 19.58 for (1, 2), N(H I) = 20.13 for (3, 
4, 5, 6, 7) and N(H I) = 19.72 for (8, 9). 
c [ergs cm" 2 s" 1 Hz" 1 sr" 1 ] 

d Density derived from the H2 as explained in Appendix 3. 



TABLE 12 

Comparison of CI ° and C II* Technique Models for DLA 1331+17 



model 


Jjocal /1Q-19 


T 


n 


log(P/k) log(%^) 




[ergs cm" 2 s -1 Hz- 1 sr -1 ] 


[K] 


[cm- 3 ] 


[cm- 3 K ] 



CII*, P b eq , min c 3.09 

CI, P e9 , min 3.09 

CII*, Pmax, min 1.74 

CI, Pram, min 1.74 

CII*, Pmin, mm 4.90 

CI, P min , min 4.90 

CII*, Peg, max d 1.64 

CI, P e q, max 1.64 

CII*, Pmax, max 0.65 

CI, Pmax, max 0.65 

CII*, Pmin, max 7.34 

CI, P m j n , max 7.34 



80.6 


6.9 


2.75 


3.70 


16-20 


229-457 


3.66-3.86 


2.40 


22.9 


17.4 


2.6 


3.23 


20-32 


95-263 


3.38-3.72 


2.40 


761.5 


1.0 


2.88 


4.51 


13-16 


347-1445 


3.74-4.26 


2.40 


102.4 


3.2 


2.52 


3.69 


25-63 


42-158 


3.22-3.60 


2.60 


23.6 


9.1 


2.33 


3.04 


25-63 


46-158 


3.26-3.60 


2.60 


838.8 


1.1 


2.97 


4.62 


13-16 


263-1096 


3.62-4.14 


2.60 



a C I results are all 2 a 
b P eq = ( P mi „ P max )V2 
c Minimal Depiction Model 
d Maximal Depletion Model 



